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Preface 


The  need  to  find  a  method  of  predicting  the  flying  qualities  of 
highly  control-augmented  aircraft  inspired  my  study  of  the  optimal  pilot 
model,  (despite  a  lack  of  knowledge  of  optimal  control  theory  and  computer 
programming) .  The  major  portion  of  the  work  done  on  this  thesis  was 
concerned  with  researching  optimal  control  theory  and  optimal  pilot 
modeling,  and  in  constructing  the  computer  program.  Difficulties  in 
obtaining  reliable  computer  service  and  my  reassignment  prior  to  the 
completion  of  this  study  limited  the  amount  of  analysis  that  could 
be  done  with  the  completed  optimal  pilot  model.  Still,  a  foundation 
has  been  established  for  further  investigation  of  this  topic  and  areas 
for  continued  study  have  been  identified. 

My  thanks  and  appreciation  are  due  to  my  thesis  advisor.  Prof. 

(Capt)  James  T.  Silverthorn,  for  his  expert  advice,  unending  patience, 
and  infinite  tolerance  of  my  slow,  rambling  progress  throughout  this  effort. 
I  would  also  like  to  thank  Capt  Richard  Floyd  for  his  help  with  the  use 
of  the  Kleinman  subroutines. 
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Abstract 

A  computer  simulation  was  constructed  using  optimal  control  theory 
to  model  the  human  pilot.  The  purpose  of  this  study  was  to  determine  if 
this  optimal  pilot  model  could  be  used  to  evaluate  highly  control-augmented 
aircraft  and  predict  the  human  pilot  opinion  rating  of  its  flying  qualities. 
Eight  aircraft/control  system  configurations  were  evaluated  twice  by  the 
model  using  two  different  sets  of  initial  conditions.  The  performance 
indices  for  each  of  the  configurations  were  plotted  and  compared  to  a 
plot  of  the  human  pilot  ratings  of  the  same  configurations .  Variances 
were  noted  in  the  relative  positioning  of  the  performance  index  ratings 
and  the  relative  positioning  of  the  pilot  opinion  ratings,  and  the  mean 
difference  in  ratings  between  the  performance  index  predictions  and  the 
pilot  ratings  varied  with  the  initial  conditions  of  the  model  (  the 
initial  conditions  of  the  actual  flight  tests  were  unknown).  The  RMS 
error  in  the  state  response  to  system  noise  for  each  configuration  was 
also  plotted  and  compared  to  the  pilot  ratings.  The  mean  difference  in 
ratings  between  the  RMS  error  predictions,  (which  were  independent  of 
the  model  initial  conditions),  and  the  pilot  ratings  compared  favorably 
with  the  mean  difference  in  ratings  between  two  human  pilots,  though  there 
appeared  to  be  no  correlation  between  the  RMS  error  and  the  pilots' 
comments  about  the  aircrafts'  flying  characteristics.  Due  to  the  small 
sample  size  involved  and  the  fact  that  the  actual  flight  evaluations  were 
not  designed  to  compliment  the  computer  simulations,  definite  conclusions 
concerning  the  viability  of  using  an  optimal  pilot  model  to  predict  pilot 
opinion  ratings  of  highly  control-augmented  aircraft  flying  qualities  could 
not  be  made;  however,  areas  for  further  investigations  were  highlighted. 
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PREDICTING  PILOT  OPINION  RATINGS  OF  FLYING  QUALITIES  OF 
HIGHLY  CONTROL-AUGMENTED  AIRCRAFT  USING  All 
OPTIMAL  PILOT  MODEL 


I  Introduction 


Background 

In  the  design  of  military  aircraft,  engineering  and  flight  ex¬ 
perience  have  shovn  that,  for  a  particular  phase  of  flight,  stability 
and  control  dynamic  characteristics  should  fall  within  certain 
specified  limits  if  an  aircraft  is  to  receive  a  favorable  pilot  rating 
of  its  flying  qualities.  These  limits  are  contained  in  the  Military 
Specification  -  Flying  Qualities  of  Piloted  Aircraft,  MIL-F-8785B, 
where  dynamic  criteria  are  listed  for  various  levels  of  flying  quality. 
Aircraft  types  and  the  phases  of  flight  are  categorized  so  that  for  a 
particular  class  of  aircraft  performing  a  specific  phase  of  flight,  the 
minimum  and  maximum  values  of  pertinent  dynamic  characteristics  can  be 
found  for  each  flying  quality  level.  Thus,  if  an  aircraft  s  designed 
to  satisfy  the  criteria  of  the  highest  level  of  flying  quality,  one 
could  anticipate  favorable  pilot  opinion  ratings  during  flight  tests. 
While  this  was  generally  true  in  the  past,  recent  increases  in  the 
capabilities  of  fighter  aircraft  and  in  the  complexity  of  their  control 
systems  have  shown  that  these  criteria  are  no  longer  always  adequate  to 
predict  flying  qualities. 

Problem 

A  study  of  highly  control-augmented  (HCA)  aircraft  was  conducted 
by  Rogers  E.  Smith  and  the  CALSPAI!  Advanced  Technology  Center  in 
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1978  (Ref.  l).  As  indicated  in  the  report,  the  YF-17  prototype  had 
been  previously  simulated  in  the  CALSPAN  ground  simulator.  The  air¬ 
craft,  an  HCA  system  which  had  been  designed  to  meet  the  criteria  in 
MIL-F-8785B,  demonstrated  satisfactory  performance  in  the  simulator. 
Actual  flight  tests,  however,  revealed  unsatisfactory  behavior  during 
the  landing  phase  and  the  pilot's  rating  of  the  YF-17  for  this  task 
was  the  worst  possible  on  the  Cooper-Harper  Pilot  Opinion  Rating 
Scale  (Fig.  l).  The  USAF / CALSPAN  variable  stability  NT-33A  aircraft 
was  flown  using  several  HCA  systems  to  gather  data  on  their  longitudinal 
flying  qualities.  Considerable  discrepancies  were  noted  when  pilot 
ratings  of  the  landing  approach  task  were  compared  to  the  flying  quality 
levels  predicted  by  MIL-F-8785B.  The  report  concluded  that  the 
specifications,  being  based  solely  on  classical  aircraft  characteristics, 
were  not  applicable  to  aircraft  with  highly  augmented  control  systems. 
Predicting  the  flying  qualities  of  these  aircraft  poses  a  problem.  A 
method  for  anticipating  aircraft  flying  qualities  that  is  applicable  to 
HCA  aircraft  is  needed. 

Objectives 

The  objectives  of  this  thesis  were  to  construct  an  optimal  pilot 
computer  program  and  evaluate  the  CALSPAN  aircraft /control  system  con¬ 
figurations  performing  an  Instrument  Landing  System  (ILS)  approach.  A 
correlation  between  the  resulting  performance  index  from  the  computer 
simulation  and  the  pilots'  ratings  of  the  actual  flight  tests  would 
indicate  the  feasibility  of  predicting  the  flying  qualities  of  HCA 
aircraft  using  optimal  control  theory. 
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Approach 

This  thesis  is  presented  in  seven  parts.  First,  Chapter  I 
provides  a  background  on  the  prediction  of  flying  qualities  and 
introduces  the  problem  of  applying  the  prediction  method  to  KCA  air¬ 
craft.  Chapter  IT  discusses  the  optimal  control  system  and  explains 
how  it  is  expanded  to  model  the  human  pilot.  Chapter  III  details  the 
modeling  of  the  aircraft  system.  In  Chapter  IV,  the  details  of  the 
construction  of  the  optimal  pilot  computer  program,  OPSACT,  is  discussed. 
Chapter  V  presents  the  results  of  the  evaluation  of  the  CALSPAK  aircraft/ 
control  system  configurations  by  OPSACT,  and  Chapter  VI  contains  the 
conclusions  and  recommendations  of  this  thesis.  The  Appendices  contain 
detailed  mathematical  developments  of  optimal  control  theory  and  the 
optimal  pilot,  an  OPSACT  program  listing,  and  the  CALSPAN  aircraft/control 
system  configuration  data. 
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x(t)  =  f3[x(t),u(t),t]  ;  x(tQ)  =  Xq  (2) 

The  performance  index  is  an  arbitrary  mathematical  expression  designed 
to  measure  how  well  a  system  is  performing  a  particular  task.  For  an 
aircraft  making  an  XLS  approach,  where  both  positive  and  negative  devia¬ 
tions  from  the  glide  slope  are  assumed  equally  undesirable,  measurement 
of  the  mean  square  error  of  the  states  and  the  inputs  from  their  nominal 
values  is  an  accurate  means  of  indicating  how  well  the  aircraft  is 
following  the  glide  slope.  Thus,  for  the  flight  task  this  thesis  inves¬ 
tigated,  the  performance  index  was 

J(u)  =  h  xT(tf)  H  x(tf) 
rx  f 

+  h  J  [x  (t)Qx(t)  +  uT(t)Ru(t) ]dt  (3) 

t 

o 

The  square  matrices  H  and  Q  are  symmetric,  positive  semi-definite 
weighting  matrices  and  R  is  a  symmetric,  positive  definite  weighting  matrix. 
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These  weighting  matrices  express  the  relative  i ’•:{«.  rt  unc-  >  :--vi  v  :  -;r.c 

cf  each  state  variable  and  control  input.  They  may  U-  exj  :  a.- 

functions  of  time.  For  example,  as  an  aircraft  near,  the  gr. tur ; 
an  approach,  deviations  in  altitude  become  more  important  re’,  r  '  . t 
other  variables.  The  element  of  the  vei.--1  '  ; matrix  Q  corre  •  i i r„-  *  c. 
the  altitude  perturbation  state  variable  shoul  1  the>’ef:  v! 4  *.  * 

Increased  weighting  indicates  to  the  optimal  control  system  4  :.  a  larger 
"penalty"  is  associated  with  deviations  from  the  nominal  alti  uue.  .Ir 
order  to  minimize  the  value  of  the  performance  index,  the  cc:  rol  system 
would  place  more  emphasis  on  controlling  deviations  in  all:-’  :e  than  on 
other  states  or  inputs.  However,  this  does  not  mean  that  .  •  '  ie  per¬ 
turbations  or  any  other  state  will  be  at  its  nominal  valut  at  4:.t.  final 
time.  The  first  term  on  the  right-hand  side  of  Eq.  ?  attempts  to  f<  r.u  s  me 
or  all  of  the  states  to  their  nominal  values  at  i.  he  final  t  ,  i>..  .  ihi>.  ih. 
however,  investigated  only  the  approach  task  and  did  not  a‘.  -  ’  to  model 

the  transition  from  the  ILS  glide  slope  to  the  landing.  The  aircraft's 
task  was  viewed  as  maintaining  an  infinitely  long  glide  si  :><r ,  thus  the 
weighting  matrices  were  constant  with  respect  to  time  and  the  first  '..era 
on  the  right-hand  side  of  Eq.  3  was  dropped. 

It  was  assumed  that  the  aircraft  system  could  be  rep:  ■  ..ented  by 
a  linear,  time  invariant  expression.  The  constraining  equal  ion  became 


x(t)  =  Ax(t)  +  Bu(t)  ;  x(tQ)  =  (M 
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u*(t)  =  -L(t)x(t)  =  -R_1BTK(t)x(t)  (5) 

where  K(t)  satisfied  the  equation 


—  K(t )  =  K(t)A  +  ATK(t)  +  Q  -  K(t)BR~1B'rK(t)  ; 

K(tf)  =  H  =  0  (6) 

Since  the  final  time  of  the  task  is  far  off  into  the  future  relative  to 
the  system  time  constants,  the  solution  to  Eq.  6  was  constant  with 
respect  to  time.  Thus  Eq.  6  became 

0  =  KA  +  ATK  +  Q  -  KBR_1BTK  (T) 

and  Eq.  5  became 

u*(t)  =  -Lx(t)  =  -R_1BTKx(t)  (8) 


Feedback  of  u*(t)  resulted  in  a  linear,  deterministic  optimal  control 
system  (Fig.  2). 


Figure  2  Optimal  Control  System 


GAE/AA/8CD-3 


Optimal  Pilot  Model 

In  1970,  Anderson  stated  that,  fiver,  a  particular  aircraft  and 
task,  the  "optimal"  pilot  will  adapt  a  control  stategv  that  minimises 
his  numerical  opinion  rating  (Ref.  2).  lie  suggested  that  combining  the 
optimal  pilot  control  strategy  with  modern  control  theory  would  yield  a 
method  of  evaluating  aircraft  performance  that  would  take  into  account 
the  pilot  and  his  rating  of  the  aircraft  dynamics.  This  combination  is 
reflected  in  the  optimal  pilot  model.  This  model  evolves  from,  the 
optimal  control  system  through  the  inclusion  of  features  that  model 
human  response  characteristics  and  the  effects  of  random  errors. 

Human  Response  Characteristics.  The  pilot  has  no  means  of 
directly  observing  the  system  states.  He  receives  state  information 
from  the  cockpit  instrumentation  and  his  perceptions  of  the  world  around 
him  (peripheral  vision  cues,  aircraft  accelerations).  If  all  infor¬ 
mation  is  assumed  to  come  from  the  instrumentation,  (as  was  the  case 
for  this  study) ,  then  the  pilot  must  estimate  the  states  based  on  the 
cockpit  display.  This  estimation  process  is  modeled  by  a  Kalman  filter 
that  reconstructs  the  states  from  the  information  available  from  the 
instrumentation.  A  detailed  mathematical  development  of  the  estimation 
process  is  contained  in  Appendix  B.  Inclusion  of  the  cockpit  instrumen¬ 
tation  and  the  estimation  process  is  depicted  in  Fig.  3. 

Human  response  includes  various  internal  time  delays.  These 
delays,  usually  due  to  lags  in  relaying  and  processing  visual  information 
in  the  brain,  are  modeled  by  lumping  them  into  a  single  internal  time 
delay,  T.  Research  has  shown  that  a  time  delay  of  0.2  seconds  is 
typical  (Refs.  3,^,5) •  The  pilot,  realizing  that  he  cannot  act  instan¬ 
taneously  upon  receipt  of  displayed  information,  attempts  to  predict  the 
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Figure  3  Modeling  Cockpit  Instrumentation  and 
the  Estimation  Process 


states  of  the  aircraft  and  compensate  for  his  inherent  time  delay.  This 
prediction  process  is  developed  in  Appendix  B  and  is  included  in  the 


model  as  indicated  in  Fig.  U  (Ref  6). 


Figure  U  Modeling  Time  Delay  and  Prediction  Process 


where  Q  and  H  are  now  symmetric,  positive  semi-definite  weightinr 
matrices  and  G  is  a  symmetric,  positive  definite  weighting  matrix. 
Ryr.ask:  and  Whitbeck  have  shown  that  adding  a  control  rate  term,  to  the 
performance  index  has  the  same  effect  as  including  a  first-order  lag 

of  the  form  1  in  the  aircraft  model  (Ref.  ?),  where  is  the 

+  1 

neuromuscular  time  constant.  However,  if  the  neuromuscular  lag  is 
modeled  simply  by  inserting  the  first-crder  lag,  the  input  from  the 
optimal  pilot  to  the  aircraft  is  calculated  incorrectly.  This  is  shewn 
below.  The  system  with  the  first-order  lag  and  an  optimal  feedback  cf 
the  aircraft  states  (~L.x)  is  depicted  in  Figure  5a.  mb.is  is  the  form  c 
the  solution  that  is  desired.  The  state  equations  become 

x(t)  =  Axx(t)  +  B  u*(t)  (10) 

where 


x(t)  = 

x(t) 

u(t) 

•  Al  = 

'a  ;  e  ' 

p  i-L 

,  and  B  = 

0 

1 

“ 

'rj 

The  optimal  input  is 

u*(t)  =  -Lx  ( t )  =Tu(t)  +  u(t)  (12) 

Y> 


If  the  augmented  state  equation  given  by  Fqn.  10  were  to  be  used  in 
conjunction  with  optimal  control  theory,  the  resulting  optimal 
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Buie  Aircraft 
Dynamic* 


(a)  Optimal  Pilot  Modal  (Open-Loop)  with  Plrat-Ordar 
taw  Added  to  Aircraft  System 


<t)  \  1 


Basic  Aircraft 
Dynamic* 


(b)  Optimal  Pilot  Model  (Closed  —Loop)  with  First-Order 
Lax  Added  to  Aircraft  System 


Basic  Aircraft 
Dynamics 


(c)  Optimal  Pilot  Model  (Closed-Loop)  with  "Equivalent* 
Alreraft/Lag  System 


!  T" 

_u_(  t) 

Basic  Aircraft 

i  Tn  «  *  1 

Dynamics 

(d)  Optimal  Pilot  Model  (Closed-Loop)  with  Equivalent  AlrcraftAag 
System  —  Input  Feedback  Loop  Combined  with  Integrator 


Basic  Aircraft 
Dynamics 


(*)  Optimal  Pilot  Model  (Closed-Loop)  with  First-Order  Lag 
and  Correct  Form  of  Optimal  Input 


Figure  5  Modeling  the  Neuromuscular  Lag 


Vw  a — 
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control  u*(t)  would  be  a?  shown  in  Figure  cb.  ’.’otic-  t*-n*  u*'4  1 
woul  i  include  not  only  the  aircraft  states  hu4  also  fee  It  a  cV  of 

innut  variable  u(t).  mhis  I.  feedback  tern,  would  effectively 

~  n+1 

nodi  fv  the  tine  constant  of  the  first-order  las'.  mo  avoid  this, 
the  following  modification  is  used. 

In  order  to  achieve  the  desired  fern  of  the  closed-loon 
optimal  input,  the  aircraft  /'lag  system,  is  modeled  as  if  u(t’*  was 
the  input  from  the  optimal  pilot,  "’his  is  done  by  rearranging  the 
expression  for  the  open-loop  input  (Fan.  12),  which  yields 


u(t) 


1  u  ( t )  =  - 1  x_(  t  ) 

^  rn 


1  u(t) 


(13) 


The  remodeled  system  is  depicted  in  Figure  he .  It  should  be  noted  that 
the  system  in  Fig.  5a  is  equivalent,  to  the  system  in  Fig.  So.  The 
open-loop  state  equation  now  becomes 


x(t)  =  A  x(t)  +  B  u(t) 
—  o~  0~ 


(1U) 


where  A  = 

r  A 

b" 

and  B  = 

’o  ’ 

o 

-  -r-- 

o 

— 

.0 

o. 

1 

L. 

do 


1? 


The  closed-loop  gain  matrix  is 
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which  is  calculated  from 


u  (t)  =  -L  x C t )  ,  L  =  P  h'h 
”  0“  c  o  o  o 


where 


Q  =  K  A  +  A  K  +Q-KB  Ff^v 
oo  oo  6  ooooo 


and  where  Q  =  T Q  !  ,  R  =  0  (where  0  is  defined  in  Fon.  Q)  . 

o  .  „  '  _  o 

.n  R. 


^he  input  feedback  loop  and  the  integrator  in  Fig.  5c  can  be  combined 

to  yield  the  system  depicted  in  Fig.  5d,  which  is  equivalent  to  the 

system  in  Fig.  5e.  From  this  figure  it  can  be  seen  that  the  correct 

feedback  gain  for  the  optimal  pilot  model  with  a  first-order  lag  is 

calculated  using  Eqns.  lU  through  l8  and  then  dropping  the  n+i,  i  =  l,2,...m, 

terms  from  L  .  ^hus,  the  optimal  feedback  gain  matrix  for  the  aircraft/ 
o 

control  system  with  the  neuromuscular  lag  is 


L  =  first  n  terms  of  L  =  I, 

1  o  - 


Therefore,  using  Fqns.  ll  -  18  to  determine  the  feedback  gain  produces  a 

closed-loop  system  as  shown  in  Fig.  Re.  ^he  desired  value  of  T  is  obtained 

n 

by  adjusting  the  values  of  the  weighting  matrices.  For  a  given 
submatrix  R  (  which  is  now  part  of  the  state  weighting 


(4 


A 
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matrix  Q  ,  where  Q  = 
o  o 

Q  !  o" 

1 

o  ;  r 

there  is  a  one-to-one  relationship  between 


R  -  G  and  T  (Ref.  6).  T  is  determined  through  an  iterative 
o  n  n 

process  of  guessing  a  value  for  R  ,  calculating  the  gain  matrix  L  , 

o  o 

comparing  the  reciprocal  of  the  n+i;  i  =  l,2,...ir.  terms  of  L  with  the 

o 

desired  values  of  ,  and  then  reguessing  Rq.  If  there  is  more  than 
one  input,  u*(t)  and  u*(t)  for  example,  the  optimal  gain  matrix  will  be 
of  the  form 


*1,1 

*1,2 

*l,n 

1 *1  ,n+l 

1 

*l,n+2 

(1?) 

*2,1 

*2,2 

*2,n 

•  *2,n+l 
i 

*2,n+2 

where  ^  n+^  corresponds  to  u*(t)  and  *2  n+2  corresponds  to  u*(t).  In 
researching  this  subject,  no  information  could  be  found  concerning  the 
consequences  of  having  more  than  one  input.  The  significance  of  the 
optimal  gains  £  5  and  i  and  the  results  of  varying  the  off- 

diagonal  terms  of  Rq  was  explored  briefly  in  this  thesis  and  will  be 
discussed  in  Chapter  IV.  More  importantly,  the  significance  of  the 
weighting  term  in  the  R  submatrix  for  the  new  state  u*(t)  and  its 
effects  on  the  optimal  gain  calculations  for  all  other  new  states 
u?(t);  i  =  2,3,...m,  was  found  to  be  critical  in  this  thesis  and  will 
also  be  discussed  later. 

Human  limitations  also  include  an  effect  called  "remnant",  (Ref.  10), 
defined  in  this  thesis  to  be  the  inherent  random  human  errors  in  the 
observation  of  the  instruments  and  in  the  execution  of  desired  control 
inputs.  The  final  human  limitation  modeled  in  this  thesis  is  "threshold". 
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which  refers  to  the  pilot's  inability  to  detect,  or  his  indifference  to, 
small  changes  in  the  instrument  gauges.  The  modeling  of  remnant  and 
threshold  is  combined  with  the  modeling  of  the  effects  of  random  external 
disturbances,  which  is  discussed  next. 

Random  Input  Effects.  Random  inputs  to  the  system  are  due  to  ex¬ 
ternal  disturbances,  such  as  wind  gusts,  and  pilot  errors.  In  investigating 
longitudinal  handling  qualities  during  an  ILS  landing  task,  only  vertical 
wind  gusts  are  considered  to  have  a  significant  effect.  Including  these 
vertical  wind  gusts  in  the  optimal  pilot  model  will  alter  the  aircraft/ 
control  system  dynamic  equations.  Eq.  10  becomes 


x(t) 


*  A  x(t)  +  B  u*(t)  +  D  w  (t) 

1—  1—  1  g 


(20) 


where  = 


D 

0 


and  w  (t)  represents  random  vertical  wind  gusts. 
8 


The  D  matrix  is  determined  by  realizing  that  the  state  aw  in  the  system 

equations  of  motion,  which  are  discussed  in  more  detail  in  Chapter  III,  is 

comprised  of  the  downward  velocity  with  respect  to  the  ground ^Wj minus 

the  movement  of  the  air  mass  w  .  For  example,  the  equation  of  motion 

8 

describing  perturbations  in  forward  velocity  is  sometimes  written  as 


Au  =  . . .  +  X  (aw)  +  . . . 

V 


Actually,  this  equation  should  be 


Au  =  ...  +  X  ( Aw-w  )+...  =  .  ..+X  Aw  -  Xv  +  ... 
w  g  w  w  g 
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Figure  6  demonstrates  the  use  of  white  noise  and  a  shaping  filter  to 
approximate  a  random  effect.  The  shaping  filter  for  verticle  wind  gusts 
at  the  altitudes  for  an  ILS  landing  approach  is  specified  in  the  military- 
specification  MIL-F-8785B  as 


(21) 


where  h^  is  the  nominal  ILS  altitude  and  Uq  is  the  nominal  forward  velocity. 
For  this  thesis,  the  wind  gust  is  approximated  for  an  altitude  of  1900  ft 
and  a  velocity  of  205  ft/sec.  MIL-F-8785B  specifies  <r  at  these  conditions 

V 

to  be  6  ft/sec.  In  this  study,  the  vertical  wind  gusts  were  modeled  as 
white  noise.  This  approximation  was  used  to  avoid  introducing  additional 
state  variables  to  account  for  shaping  filters.  This  avoided  the  compu¬ 
tational  complexities  of  a  large  number  of  variables.  From  Ea.  21,  the 
power  spectral  density  of  the  vertical  wind  gusts  is  flat  with  an  amplitude 
given  by 


T  (s)  =  18.3 

g  s=0 


Thus,  the  intensity  of  the  gusts  is 


W  =  (18. 3)2  =  33U 


Human  controller  remnant,  as  defined  by  Levison,  Baron,  and  Kleinman  (Ref.  10), 
is  the  portion  of  the  pilot's  input  not  associated  with  the  aircraft/ 
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control  system  input.  Though  vertical  wind  gusts  are  normally  included 
as  part  of  the  system  dynamic  equations,  this  thesis  models  them  along 
with  observation  and  motor  noise  as  a  separate  input  source. 

Motor  noise  consists  of  the  neuromuscular  errors  the  pilot  makes 
in  activating  the  controls  and  accounts  for  the  fact  that  he  does  not 
have  perfect  knowledge  of  his  control  inputs.  In  modeling  the  effects 
of  motor  noise,  it  should  he  recalled  that  though  the  aircraft /control 
system  had  to  he  modeled  according  to  Eq.  14  in  order  to  generate  the 
appropriate  feedback  gains,  the  system  still  has  the  same  control 
structure  and  is  modeled  according  to  Eq.  10  (Ref.  6).  By  including  the 
motor  noise,  the  control  input  defined  by  Eq.  8  becomes 

u*(t)  +  vu(t)  =  -Lx(t)  +  vu(t)  =  Tu(t)  +  u(t)  (22) 


or  equivalently 


u(t) 


fn  x(t)  + 


V  (t) 
u 


\J  {  i  ! 

T 

'  -I TV 


(23) 


where  v  (t)  represents  the  motor  noise.  The  system  is  now  modeled  by 
including  the  motor  noise  into  Eq.  20 


“x(t)  =  A-.'x(t)  +  B  u*(t )  +  D  v  (t) 

X  ±  C.  g. 


(2k) 


where  v„  (t)  = 
gl 


vg(t) 


v  (t) 
u 


and  = 


D1  i  0 

o  !  i 
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-r- 


GAE/AA/80D-3 


Studies  of  human  remnant  have  shown  that  motor  noise  is  pro¬ 
portional  to  the  variance  of  the  inputs  (Ref.  10).  Specifically, 


V  =  .003U'  E 
u 


■|  [u*(t)-u* 


,*(t)][u*(t)-u*(t)] 


■} 


(25) 


Also,  the  intensity  of  the  external  disturhance/motor  noise  is 


Wl  =  D2E{VtV<t)} 


T  1 

D  WD  I  0 
1  1  I 


- 1 - 

o  ,  vu 

1  r  2 

I  'n 


(26) 


The  determination  of  the  motor  noise  matrix  V  is  done  in  an  iterative 

u 

manner.  A  guess  for  must  be  made  in  order  to  form  W^.  The  estimator 
gain  is  then  calculated  (Appendix  B).  The  steady  state  covariance  of  the 
state  vector  x(t),  which  contains  the  variance  of  the  subvector  u*(t),  is 
then  calculated  using  a  relationship  developed  by  Kleinman  (Ref.  ll). 

This  relation  is 


{rrl  A,T  A,T  /•-  ft  (T  A  U 

(x-^)(x-x)Ti  =  e1Pe1  +  J  {e  -  W  e  1  ]  d«T 
'  o 

rm  A<r  at  Air  a«t 

+  J  [e  e  PC  V  CP  e  e  ]d<r 


-T  A1<r  aV 


(27) 


where  A  =  A  -B  L  ,  P  is  the  estimator  error  covariance  matrix  (to  be 
o  o  o 

discussed  later),  and  is  the  observation  noise  matrix  (to  be  discussed 
next).  The  n"^  +  1  element  on  the  principle  diagonal  of  the  X  matrix. 
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(corresponding  to  the  n  “'  +  1  state,  u*(t)  ;  i  =  l),  is  the  variance  of 

2  2 
the  input,  <P  .  The  first  guess  of  V  is  now  compared  to  .003*®'  * 

U  Ui,i  Un+i 

If  the  difference  is  greater  than  some  arbitrary  tolerance,  (5%  of 
2 

.003ir<T  was  used  for  this  thesis),  then  the  new  guess  of  V  is  set 
u  .  .  &  u.  . 

n+i  2  1,1 

equal  to  .003*®  and  the  procedure  is  repeated.  This  process  is 

Un+i 

described  in  more  detail  in  Chapter  IV. 

The  final  component  of  remnant  modeled  in  this  thesis  is  the 
observation  noise.  This  effect  accounts  for  the  errors  the  pilot  makes 
in  reading  the  instruments  and  the  random  errors  in  the  instruments 
themselves.  Similar  to  the  motor  noise,  the  observation  noise  is 
proportional  to  the  variance  of  the  displayed  variables.  Specifically , 

Vy  =  .01*  E|[£(t)-jr(t)][y(t)~£(t)]Tj  (28) 


Determination  of  V  is  done  in  the  same  iterative  manner  as  with  V  , 

y  u 

noting  that 


E  {[jl(t)-JL(t)][i(t)-i(t)]T  | 


|[x(t)-x( 


=  C  E  i  [x(t)-x(t)][x(t)-x(t)]T |  CT  (29) 


where  C  is  the  displayed  variable  (instrumentation)  coefficient  matrix. 

One  other  item  is  included  in  the  determination  of  the  obser¬ 
vation  noise  matrix.  When  reading  the  instrument  gauges,  the  pilot 
cannot  detect  changes  in  the  readings  that  are  smaller  than  a  certain 
threshold  value.  Also,  the  pilot  may  choose  to  ignore  readings  that  he 
can  see  but  are  smaller  than  his  indifference  threshold.  Since  the 
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I 


♦ 


♦ 


indifference  threshold,  differs  from  pilot  tc  pilot,  it  was  assumed  that 
the  optimal  pilot  would  react  to  any  changes  ir.  the  gauges  that  he  could 
detect.  The  visual  threshold  was  modeled  assuming  that  the  pilot  looks 
directly  at  the  instrument  panel.  Studies  have  shown  that  when  viewing 
an  object  straight  on,  the  human  eye  can  detect  changes  in  position 
equivalent  to  .05°  of  visual  arc  and  changes  in  rates  equivalent  to  .18° 
of  arc  per  second  (Ref.  3).  ''or  example,  if  a  vertical  velocity  indicator 
was  3  inches  in  diameter  and  was  calibrated  as  shown  in  Pig.  7,  a  movement 
of  the  indicator  of  .02b  inches  would  correspond  to  1  ft/sec  in  vertical 


Figure  7  Determination  of  Visual  Threshold 


i 
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Assuming  the  instrument  was  30  inches  away  from  the  pilot,  the  visual 
arc  swept  out  by  the  .02^  inch  movement  of  the  indicator  is  .0^5°/ft/sec 
Dividing  the  minimum  visual  arc  perceivable  by  the  actual  arc  swept  ou* 
results  in  the  visual  threshold  J\  In  this  example,  3*=  1.11  ft/sec. 
Therefore,  any  vertical  velocity  indication  less  than  1.11  ft/sec  would 
not  be  detected  by  the  pilot. 

The  information  the  pilot  is  able  to  perceive  y  (t)  is  now 

P1 

described  by  the  non-linear  function  f(y)  and  the  observation  noise 
vy(t). 


y  (t)  =  f(y)  +  v  (t) 

Ft  J 


(30 


The  function  f(y)  is  depicted  in  Fig.  8. 
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In  order  to  incorporate  the  threshold  effect  into  the  pilot  model,  a 
linear  representation  of  f(y)  is  needed.  This  representation  is 
developed  by  Kleinman  (hef.  3).  A  difference  function  d(t)  is 
defined  as  follows: 


d(t)  =  f(y)  -  f  •  y(t) 


where  f-y  is  the  linear  approximation  of  f(y)  and  minimizes  the  relation 
T  =  E  jd2(t)j  =  E  jjjf ( y )  ~  f  •  y(t^| 

■  E|[f(y)]2  -  2f  *  y ( t ) [ f ( y ) ]  +  f2  •  y2(t)}  (32) 


To  minimize  T,  Eq.  33  must  be  true. 


0  =  E  [0]  -  2y(t)f(y)  + 


A 

2[f 


y2(t )] 


(33) 


Thus , 


A 

f 


=  E  <  y(t)f(y) 


y(t)f(y)p(y)dy 


2 

y 


(t)p(y)dy 


-1 


(3U) 


Since  y(t)  is  assumed  to  be  a  Gaussian  random  variable. 


p(y)  =  probability  density  function  of  y 
1  1 


-(y-  y)2/2C2 


e 


(36) 
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where  cr=  the  standard  deviation  ar.  y  =  the  near..  Substituting  p(y) 

A. 

into  Eq.  3t  and  assuming  y(t)  is  a  zero-mean  process,  f  becomes 


S 


(36) 


where  er 
becomes 


fc(*y2} 


is  the  complementary  error  function  of 


Eq.  30 


yp  (t)  -  yp  (t)  =  f  •  y(t)  +  Vy(t) 


(37) 


Since  the  pilot  will  mentally  correct  for  this  threshold  effect,  the 
measure  he  uses  is 


yp  (t)  =  yp  (t)  =  y(t)  +  vy(t'  (38) 

3  _ 2 _  , 


Since  0<f— 1,  -this  has  the  effect  of  increasing  the  noise  present  in  the 

A_1 

measurement  by  a  factor  f  .  Therefore,  to  correct  for  this,  Eq .  26  becomes 
Vy  =  -Olir  E  |[x(t)-^(t)][^(t)-£(t)]T  }  •  f'2  (39) 


Problems  arise  in  calculating  the  threshold  effects  if  the  chosen  visual 

A 

threshold  level  f  is  much  larger  than  <Ty.  As  can  be  seen  in  Fig.  9,  f  is 
represented  by  the  shaded  area  under  the  Gaussian  distribution  curve,  (the 

A 

total  area  tinder  the  curve  is  equal  to  1.0).  If  £>>0^,  then  f  becomes 

A-1 

very  small,  and  f  gets  very  large.  Figure  10  shows  the  relationship 

A  A 

between  f ,  #,  and  <Ty.  If  the  value  of  df  is  such  that  <Ty/(f  •  f  )  -*  00 
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then  the  solution  to  the  iterative  V  determination  process  will  not 
converge.  Hence,  the  proper  calculation  of  the  visual  threshold  value 
is  critical. 


p(y) 


Figure  9  Complimentary  Error  Function 


Figure  10  The  Effects  of 


Summary .  At  this  point,  the  portion  of  the  optimal  pilot  model 


that  deals  strictly  with  the  human  pilot  has  been  developed  from  the 
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optimal  control  system.  Human  response  characteristics,  (state 
estimation,  internal  time  delay,  neuromuscular  lag,  visual  threshold, 
remnant),  and  the  effects  of  random  external  disturbances  have  been 
modeled.  The  proper  methods  of  determining  the  optimal  feedback 
gain,  observation  and  motor  noise  matrices  have  been  discussed. 

What  remains  of  the  optimal  pilot  model  to  be  explained  is  the  modeling 
of  the  aircraft /control  system.  This  explanation  follows  in  Chapter  III. 


GAE/AA/80D-3 


III  Modeling  the  Aircraft/Ccntrol  System 


Aircraft/Control  System 

The  figure  below  depicts  the  major  components  of  the  aircraft  system 
The  discussion  of  the  modeling  of  the  aircraft/control  system  will  address 


each  of  these  components. 


Figure  11  The  Aircraft  System 


Basic  Aircraft  Dynamics.  The  basic  aircraft  equations  of  motion 
are  (Ref.  l) 
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It  should  be  noted  that  the  equations  of  motion  given  by  Ref.  1  have  been 
rearranged  in  order  to  imbed  the  M.  and  Z.  terms  within  the  remaining  terms 
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Thus,  the  stability  derivatives  above  are  actually  equivalent  values  and 
include  the  effects  of  and  Z^.  The  pilot  is  unable  to  directly  observe 
most  of  the  variables  listed  above.  Though  the  cockpit  instrumentation 
will  display  some  of  these  variables,  most  of  the  pilot  information  about 
the  aircraft  will  consist  of  linear  combinations  of  the  actual  states. 

One  important  item  of  information  is  the  cockpit  display  of  the  error  in 
tracking  along  the  ILS  glide  slope.  The  error  is  in  the  form  of  At  ,  the 
angle  above  the  nominal  glide  path  angle  f^.  In  order  to  relate  AC  to 
the  dynamics  of  the  aircraft,  it  is  necessary  to  add  additional  states 
to  the  basic  aircraft  dynamic  equations.  These  states  are  perturbations 
in  altitude  Ah  and  perturbations  in  slant  range  Ar.  The  states  and  coordinate 
systems  are  depicted  in  Fig.  12.  From  the  figure  the  following  can  be  seen: 


A 

e 


3 


(-sin  6)  +  (cos  e) 


(in) 


e  „  =  -cos 

gs 


(T 


e)  ^  -  sin  (T0  +  e) 


A 


(h2) 


Thus, 


h  =  U  sin  0  -  W  cos  6  =  h^  +  Ah  =  f^fU,  W,  6) 

r  =  -U  cos  (P  +  0)  -  W  sin  (T  +  S) 
o  o 

=  r  +  Af  =  f  (u,  w,  e) 

O  d 


<U3) 


(UU) 


The  perturbation  equations  are 
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Ah 


£ 

bu 


Au  + 


tv 


Aw  + 


*fl 

be 


£»© 


=  Au  sin  0Q  -  Av  cos  0^ 


+  (U  cos  0  +  W  sin  0  )Ae 

o  o  o  o 


(1*5) 


af2 

bf2 

. 

Ar  =  - 

Au  +  — - 

Aw  +  — - 

au 

aw 

o 

Ae 

o 

A© 


=  -A u  cos  ( T  +  6  )  -  aw  sin  ( T  +  6  ) 
o  o  o  o 


+  a©[Uq  Sin  (  T0  +  Bq)  -  VQ  cos  (  Pq  +  6Q)]  (1*6) 


It  is  reasonable  to  assume  that  the  nose  of  the  aircraft  is  pointed 
approximately  along  the  glide  slope;  thus, 


e„  =  -r 


0*7) 


Equations  1*5  and  1*6  become 


Ah  =  -Au  sin  P  -  Aw  cos  P 
o  o 

+  (U  cos  P  -  V  sin  P  ) A© 
O  0  0  o 


(1*8) 


A  r  =  -Au  -  V  A© 

o 


(1*0) 


The  nominal  glide  path  angle  used  was  .01*5  rads  (2.5°);  thus,  was 
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assumed  to  be  a  small  angle.  Equations  ^8  and  U9  become 


Ah  = 


-auP  -  Aw  +  (U  -  W  f  )Ae 

o  o  o  o 


Ar  =  -Au  -  W  Ae 

o 


The  basic  aircraft  dynamics  become 
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Stick  Feel  System/Actuator.  The  dynamics  for  the  control  stick 
feel  system  and  the  actuator  for  the  NT-33A  are 


Se. 


8^.5 


Feel  system: 


*e. 


s  +  36. Us  +  676 
5625 


Actuator: 


cs 


+  105s  +  5625 


These  systems  were  dropped  from  the  model  because  of  their  high  response 
frequencies.  It  was  assumed  that  for  the  human  controller  problem  the 
frequencies  of  interest  would  be  in  the  0  to  l6  rad/sec  range. 
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Control  System.  Various  control  systems  (listed  in  Appendix  D) 
were  tested  in  this  thesis.  To  demonstrate  how  the  basic  aircraft 
dynamics  were  augmented  by  these  control  systems,  a  typical  configuration 
is  chosen  as  an  example.  Control  System  7  is  described  by 

gecs  i1*1* 

&ec  s2  +  16.8s  +  l4U 

In  phase  variable  form  this  becomes 

X1  =  ^c*  ;  X1  =  x2 

* 

Xg  =  ;  x2  =  -l6.8  x2  -  lUU  x.^  +  lUU 

The  aircraft  dynamics,  including  the  effects  of  external  disturbances  and 
motor  noise  as  indicated  by  Eq.  2h,  are  now  described  by  Eq.  53  on  the  next 
page. 

An  additional  control  input,  thrust  control  ,  was  added  because 
the  elevator  alone  will  not  adequately  control  the  aircraft  on  the  glide 
slope.  This  throttle  term  comes  from  the  realization  that  the  thrust  forces 
acting  on  the  aircraft  are 

=  Th  cos  l/  +  Th  siny 

where  V  is  the  angle  between  the  engine  thrust  axis  and  the  aircraft 
centerline.  Assuming  that  V  =  0  and  that  there  are  no  lateral  forces  due 
to  thrust,  the  forces  acting  on  the  aircraft  due  to  thrust  are 
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Th 


Th 


Adding  the  perturbation  thrust  to  the  Au  equation  of  motion  yields 


Au  =  X  Au  +  X  Aw 
u  w 


W  AQ  +  X< 
o  *  *«  *es 

cs 


ATh 

m. 


X  v 
w  g 


(5M 


For  an  ideal  jet,  ATh  is  dependent  only  on  altitude  and  power  setting. 
For  the  small  altitude  changes  examined  in  this  thesis,  ATh  is  dependent 
only  on  the  power  setting  £  .  Thus, 


ATh 

“l 


(5J 


where  Th  is  the  maximum  thrust  available  and  jtl  is  the  weight  of  the 
max  1  c 

aircraft.  For  the  NT-33A,  Th  is  approximately  5 >000  lb„  and  the 

max  i 

weight  for  the  approach  is  approximately  12,500  Ib^.  Thus,  the 
coefficient  for  the  5^  term  is 


For  this  thesis,  the  state  variable  Ar  was  not  included  in  the 
modeling  of  the  aircraft /control  system.  The  optimal  pilot  computer 
program  subroutine  designed  to  solve  Eq.  7  (subroutine  MRIC  in  Appendix  C) 
would  not  converge  towards  a  solution  if  this  state  was  included  in  the 
model.  No  reason  for  this  difficulty  could  be  determined,  so  the  Ar 


state  was  dropped  and  the  value  of  r  for  all  applicable  terms  in  the 
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displayed  variable  coefficient  matric  C  (to  "be  discussed  next)  was 
approximated  as  r  ,  the  nominal  slant  range. 

Cockpit  Instrumentation.  As  previously  mentioned,  the  pilot  is 
unable  to  observe  most  of  the  aircraft  states  directly.  The  cockpit 
instrumentation  displays  some  of  the  states  directly,  but  most  of  the 
information  displayed  is  a  linear  combination  of  the  states.  The 
instrumentation  available  to  the  NT-33A  pilot  for  the  glide  slope 
task  is  depicted  in  Fig.  13-  It  is  usually  assumed  that  if  the  pilot 
can  read  a  displayed  quantity  explictly,  then  he  can  implictly  obtain 
information  about  the  rate-of-change  of  that  quantity.  However, 
instrument  lag  due  to  hysteresis  causes  doubt  concerning  the  usefulness 
of  implicit  rate  information  during  an  ILS  landing  approach.  Therefore, 
the  displayed  variables  modeled  in  this  thesis  were 

(1)  Perturbation  angle  of  attack  (&“) 

(2)  Perturbation  rate-of-change  of  angle  of  attack  (&<*) 

(3)  Perturbation  pitch  attitude  (a9) 

(U)  Perturbation  rate-of-change  of  pitch  attitude  (A9) 

(5)  Perturbation  airspeed  (&<*) 

(6)  Perturbation  descent  rate  (aM 

(7)  Perturbation  angle  above  the  nominal  glide  slope  angle  (AC} 
Those  displayed  variables  that  are  linear  combinations  of  the  system 
states  are  developed  as  follows: 

1 

A«  =  U  Aw  (5  6) 

o 

1  z  z  Z$e0 

A*  =  —  Aw  =  ~  Au  +  Aw  +  Aq  +  — - A  ( 57 ) 

o  o  o  o 

Ah  =  -T0Au  -  Aw  +  (Uq  -  W oro)Ae  (58) 
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The  display  variable  is  developed  as  follows: 


£  =  €  +  A€  ;  £  =  0 

o  o 


o  r 

o 


The  displayed  variables  are  now  determined  according  to  Eq.  60  on  the 
next  page. 

Weighting  Matrices.  Determination  of  the  values  of  the  weighting 
matrices  is  a  difficult  task.  While  the  aircraft  dynamics  can  be  calculated 
from  experimental  and  analytical  data,  the  calculation  of  the  weighting 
matrices  is  more  subjective.  Hess  suggested  that,  in  order  for  an  optimal 
pilot  model  to  predict  the  pilot's  rating,  the  elements  of  the  weighting 
matrices  should  be  the  reciprocal  squares  of  the  maximum  desirable  pertur¬ 
bations  of  their  corresponding  variables  (Ref.  7).  The  change  in  any 
variable  would  be  viewed  as  a  fraction  of  the  maximum  desirable  perturbation. 
The  penalty  assessed  by  the  performance  index  for  a  fractional  change  in 
one  variable  would  be  equal  to  the  penalty  for  the  same  fractional  change 
in  another  variable.  For  example,  consider  two  input  variables; 
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elevator  perturtat ions  and  thrust  perturbations  .  Figure  ll  shows 
that,  if  the  weightings  are  chosen  as  Hess  suggested,  a  change  of  lOf  of 
the  maximum  desirable  perturbation  of  A£e  would  ’’ncur  the  same  penalty  as 
a  1C:"-  change  in  The  optimal  controller,  in  attempting  to  minimize 

the  penalties  assessed  against  all  variables,  will  not  favor  one  variable 
over  another,  (though  it  would  if  the  penalty  for  a  fractional  change  in 
one  variable  were  more  than  for  the  same  fractional  change  in  another 
variable) . 

Difficulty  arises  in  determining  the  value  of  the  maximum 
desirable  perturbation.  Obviously,  these  values  will  differ  from  pilot 
to  pilot.  This  means  that  a  choice  of  weighting  values  would  yield  an 
optimal  pilot  model  that  would  predict  the  ratings  of  only  the  pilot 
whose  subjective  weighting  values  were  used.  It  was  assumed,  however, 
that  well-trained,  skillful  pilots  do  not  differ  significantly  in  what 
they  consider  acceptable  control  and  state  perturbations  during  an  approach, 
and  that  these  differences  would  have  very  little  effect  on  the  results  of 
the  optimal  pilot  model. 

In  choosing  the  weightings  for  this  thesis,  the  opinions  of  the 
CALSPAII  pilots  who  flew  the  actual  missions  were  not  available,  so  the 
opinions  of  available  pilots  experienced  in  fighter  aircraft  were  used. 

These  weightings  are  depicted  in  Fig.  15 . 

It  was  felt  that  the  pilot  cannot  make  reasonably  accurate  assess¬ 
ments  of  the  maximum  desirable  perturbations  of  Aq  and  ,  so  the 

weightings  were  set  to  zero.  This  would  have  the  effect  of  telling  the 
optimal  controller  that  perturbations  in  these  quantities  are  not 
important  to  the  pilot.  Actually,  by  weighting  other  variables,  these 
variables  will  remain  well-behaved.  The  weighting  on  was  originally 
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set  to  131-3  (5°  perturbation).  This  was  later  set  to  zero  to  allow  the 
desired  value  for  the  neuromuscular  lap  to  be  achieved.  This  is  discusse 
further  in  Chapter  IV. 

Summary.  The  aircraft/control  system  modeling  has  been  explained 
and  the  developement  of  the  optimal  pilot  model  is  now  complete.  Chapter 
IV  will  now  discuss  the  organization  and  construction  of  the  optimal 
pilot  computer  program,  OPSACT. 
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A  Aircraft/Control/Equivalent  Lag  (A/C/KL)  System  Coefficient 

Matrix. 

A/C/EL  System  Input  Coefficient  Matrix. 

C  Displayed  Variables  (Instrumentation)  Matrix. 

Qq  A/C/EL  System  State  Weighting  Matrix. 

Rq  A/C/EL  System  Input  Weighting  Matrix. 

V  Displayed  Variable  (Observation  Noise)  Covariance  Matrix 

W^  External  Disturbance/Motor  Noise  Covariance  Matrix 

t  Initial  Time  of  Task 

o 

t_  Final  Time  of  Task 

At  Integration  Time  Increment 

x  Initial  State  Vector 

o 

uq  Initial  Input  Vector 

Jf  Visual  Tresholds 

Figure  IT  OPSACT  Inputs 

To  solve  Eq.  10,  OPSACT  integrates  from  tQ  to  tf  using  transition 
matrices.  Eq.  10  becomes 

A  At  -(k+l)At  a  [  (k+ijAt-V] 

xJ(k+l)At]  =  e  x(kAt)  +  J  e  B  u (C ) d IT 

kAt  1 

=  <|)(At)j[(kAt)  +  F(4t)u(kAt)  (6l) 


The  estimator  equation,  (Eq.  B-6,  Appendix  B),  becomes 
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A 

x[(k+l)At]  = 


+ 


A  it* 

e  1  x(kAt) 


*  / 


(k+1 )At 


kAt 


A  [  (k+1) At -O'] 
e  B^u(<r)d<r 


f 


( k+1 ) At 


kAt 


A  [(k+1) At a 
e  S  [  y_(  kAt )  -Cx  ( k  At )  ]  d<T 


A  " 

=  $(At)x(kAt)  +  f1  (At  )u(kAt )  +  4,(At)[jr(kAt)-Cx(kAt)] 


(62) 


The  predictor  equation,  (Eq.  B-20),  becomes 


xt(k+l)At] 


A.,  At* 

e  x*(kAt)  + 


kAt 


A 

=  <()  ( At  )x* (kAt )  +  T(At)u*(kAt) 


(63) 


A  At 

Since  At  =  constant  for  this  thesis,  <t>(At)  =  e  =  constant.  By  making 
the  transformation  £=  [  (k+l)At-<7] ,  then  T(At)  and  ^  (At )  become 


f{ k+l)At  A  [(k+l)At-<T]  A J 

r(At)  =  J  e  B  d<T  =  J  e  1  djB 

kAt  o 


=  constant 


-(k+l)At  A  r(k+1)At-ff] 

At  . 

r 

f(At )  =  J  e  Sd  = 

kAt 

(6U) 


=  constant 


(65) 


A1T 

Thus,  <f>,  P,  and  another  constant  used  by  the  predictor,  e  ,  are 
calculated  outside  of  the  pilot/aircraft  integration  loop.  The  estimator 


1*3 
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gain  matrix  S  is  needed  from  part  three  of  the  program  before  the 

calculation  of  is  complete. 

Optimal  Feedback  Gain/Neuromuscular  Lag 

The  optimal  feedback  gains  are  determined  according  to  Eqs .  IT  and  18. 

The  iterative  process  of  arriving  at  a  desired  value  for  the  neuromuscular 

lag  Y  through  repeated  guesses  of  R  and  P.Q  ,  as  discussed  in  Chapter  II, 

can  be  accomplished  by  reentering  OPSACT  with  different  values  for  R  and 

Rq  .  This  method,  however,  consumes  a  great  deal  of  computer  time.  At 

present,  there  is  no  feature  in  OPSACT  to  allow  the  entering  of  desired 

values  of  T  and  letting  the  computer  vary  R  and  R  in  order  to  obtain 
n  o 

the  desired  lag.  For  this  thesis,  a  second  program,  OPTCOE,  (which  is 
designed  to  solve  only  the  optimal  feedback  gain  problem),  was  used  to 
obtain  the  correct  values  of  R  and  Rq  that  would  yield  the  desired  neuro¬ 
muscular  lag  for  each  aircraft/control  system  configuration.  The  inputs 

to  the  optimal  feedback  gain  problem  are  A  ,  B  ,  Q  and  R  .  No  problems 

o  o  o  o 

were  encountered  in  arriving  at  the  desired  thrust  neuromuscular  lag 

of  1.0  second,  (this  lag  actually  included  the  engine  lag  along  with 

the  pilot  lag) .  Considerable  difficulty  arose  in  finding  values  of 

R  and  Rq  that  would  yield  the  desired  elevator  lag  of  0.1  seconds.  Results 

showed  that  as  the  corresponding  element  in  R  decreased,  or  the  corresponding 

element  in  Rq  increased,  the  calculated  lag  would  increase.  Initial 

guesses  of  R  and  R  resulted  in  too  small  a  value  for  Y  .  Repeated  iterations 
o  n 

e 

resulted  in  the  R  element  being  reduced  to  zero  and  the  R  element  being 

o 

increased  to  the  point  that  0PTC0N  could  not  arrive  at  a  solution,  yet 

the  value  of  Y  was  still  too  small.  Experimenting  with  the  off-diagonal 
e 

elements  of  Rq  revealed  that  any  non-zero  entry  in  these  elements  would 
only  decrease  the  value  of  Y^  and  Y  .  It  was  found,  however,  that  by 
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increasing  the  weighting  on  the  commanded  throttle  input  in  the  F: 

c 

submatrix  from  25.0,  (which  corresponded  to  a  maximum  desirable  perturbation 
Tfr 

of  20$  ) ,  to  a  weighting  of  1*00.0,  (5$  ~ —  ),  and  then  reguessing 

max  max 

R  to  obtain  the  1.0  second  thrust  lag,  the  element  of  R  could  now  be 

increased  to  achieve  the  desired  elevator  lag  without  having  computational 

problems.  This  relationship  between  the  weighting  of  one  input  and  the 

arrival  at  the  desired  neuromuscular  lag  of  another  input  was  not 

mentioned  in  any  of  the  research  material. 

The  final  action  taken  by  the  second  part  of  OPSACT  is  to  take 

the  calculated  values  of  and  form  the  actual  system  state  coefficient 

matrix  A_^  as  described  by  Eq.  10.  The  values  of  T^  are  then  stored  for 

use  in  forming  the  external  disturbance/motor  noise  covariance  matrix  W  . 

Estimator  Gains/Variances 

The  estimator  gain  matrix  S  is  calculated  using  the  A^  matrix. 

The  third  part  of  OPSACT  forms  the  matrix  and  determines  the  estimator 

gains  according  to  Eqs.  B-17  and  B-l8.  As  explained  in  Chapter  II,  the 

motor  noise  V  is  proportional  to  the  variance  of  the  n+i  states  u.(t); 
i 

i  =  l,2,...m;  and  the  observation  noise  V  is  proportional  to  the  variance 

yi 

of  the  displayed  variables  y^(t);  i  =  l,2,...ny,  where  ny  =  the  number 
of  variables  displayed  on  the  instrument  panel.  The  covariance  matrix  of 
the  system  states  is  calculated  using  Eq.  27-  The 'first  two  terms  of  that 
equation  are  calculated  as  indicated.  Noting  that  the  third  term  can  be 
written  as  (Ref.  lU) 


A°*  AT  A^T  AV 

J  e  [e  1  pcV^CP  e  1  ] 
o  y 


tT_ 

e  a<r= 


Q 

f 


at  aV 

e  G  e  (66) 
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from  part  2 


to  part  4 


Figure  18  Part  Three  of  OPSACT 
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Evaluation  of  the  integral  in  Eq.  66  is  equivalent  to  finding  the 
solution  matrix  Z  to  the  equation 


0  =  AZ  +  ZA1  +  C- 


(67) 


Eq.  67  is  a  Lyapunov  equation  and  is  solved  to  determine  the  third 
term  of  Eq.  27.  The  covariance  matrix  of  the  displayed  variables  is 
obtained  from  the  results  of  Eq.  27. 


T 

Y  =  C  XC 


(68) 


As  depicted  in  Fig.  19,  the  n=i  diagonal  element  of  V  ;  i  =  l,2,...m, 

2  2  2 

is  multiplied  by  T  and  is  compared  to  .  0031KJ*  ,  where  is  the 

n.  u.  u. 

l  li 

n+i  diagonal  element  of  the  state  covariance  matrix  X.  If  the  difference 

2  2 
is  greater  than  5%  of  .003VO  ,  then  V  is  set  equal  to  .003^  and 

2  i  ^ 

divided  by  *T  to  form  a  new  matrix.  Similarly,  the  i  diagonal  term 

ni  th 

of  V^;  i  =  l,2,...ny,  is  compared  to  the  i  diagonal  term  of  Y, 

.01'itG‘^/f^.  If  the  difference  is  greater  than  5$  of  .0r*r<f*  ^/f" ,  then  V 
yi  yi  y 

2  2 

is  set  equal  to  .OlTTC^/f  .  If  either  or  was  changed,  the  program 
returns  to  the  start  of  part  three  and  repeats  the  process  until  all 
elements  of  V  and  V  are  within  5%  of  the  desired  values.  The  program 

u  y 

then  continues  to  part  four. 

Pilot/Aircraft  Integration  Loop 

Part  four  of  OPSACT  combines  the  actions  of  the  aircraft  system, 
the  estimator,  and  the  predictor  to  form  the  pilot/aircraft  iteration  loop. 
This  process  divides  the  continuous  actions  of  the  pilot/aircraft  system 
into  small  increments,  where  each  increment  represents  one  pass  through  the 


L7 
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It  should  be  noted  that  the  order  in  which  the  calculations  are  made  is 
not  determined  by  following  the  arrows  around  Fig.  20.  Though  the  looj 
models  the  control  task  in  discrete  steps,  the  actual  integration  process 
made  by  the  pilot/aircraft  system  is  continuous.  Fig.  20  only  indicates 
where  the  various  inputs  and  outputs  go.  The  model  is  not  a  sequential 
system  that  is  initiated  at  the  aircraft /control  system  block  with 
the  input/output  information  flowing  along  the  arrows  activating  other 
components  of  the  model.  During  any  one  pass  through  the  loop,  inputs 
exist  for  all  of  the  integrations,  (initial  values,  new  values,  delayed 
values);  thus,  the  order  in  which  the  integrations  take  place  is  not 
important.  Figure  21  shows  which  components  must  precede  others, 

(for  example,  the  delaying  of  x(t)  must  follow  the  determination  of 
that  quantity).  The  grouping  of  the  integration  processes  indicates 
that  they  should  be  viewed  as  taking  place  simultaneously. 

The  choice  of  the  integration  time  increment  At  is  not  totally 

arbitrary.  When  attempting  to  reconstruct  a  signal  from  discrete 

samples,  a  sampling  rate  that  is  too  slow,  (time  increment  too  large), 

will  prevent  accurate  reconstruction  of  the  original  signal.  Figure  22 

shows  that  if  a  signal  Sq  is  sampled  at  a  low  frequency  f ^ ,  the  reconstructed 

signal  S  will  not  resemble  the  original  signal.  Sampling  at  a  higher 
rl 

rate  f  ^  will  yield  a  reconstructed  signal  that  is  a  better  represen4 ation 
of  the  original  signal.  Shannon's  sampling  theorem  states  that  a  signal 
can  be  accurately  reconstructed  if  the  sampling  increment  is  chosen  such  that 


At 


h9 


T"* 


Pigure  21  Order  of  Calculation  of  Components 
of  Pilot/Aircraft  Integration  Loop 


Figure  22  Signal  Sampling  Rate 


where  is  the  highest  frequency  in  the  signal  (Ref.  12).  In  practice, 
other  factors  require  that  a  higher  sampling  rate  he  used.  Common  choices 
are  6  to  20  samplings  per  signal  cycle.  For  this  thesis,  this  would  mean 


1  21T 


that  the  integration  time  increment  would  have  to  he  between  y  —  and 

6  to 


1  2T  , 

—  —  ,  where  is  the  largest  eigenvalue,  ^highest  frequency  m  the  system 
c 

response).  Since  the  pilot's  internal  time  delay  T  is  achieved  by 


OPSACT  by  delaying  certain  quantities  an  integer  number  of  loop  iterations. 
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At  was  chosen  as 


4ts  5  ('»*)  =  .025  seconds  =  J  (0.2)  =  ^  (£?)  (70) 

c 

where  o>c  was  chosen  as  17  rads/sec,  which  was  the  highest  frequency  in 
the  response  of  the  systems  tested.  This  made  the  sampling  rate  approx¬ 
imately  15  times  per  second. 

A  listing  of  program  OPSACT  is  contained  in  Appendix  C.  Included 
in  this  listing  are  the  subroutines  written  by  Kleinman  (Ref.  lc). 
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IV  Results 


In  order  to  predict  the  flying  qualities  of  aircraft,  a  correlation 
must  exist  between  the  performance  index  generated  by  the  optimal  pilot 
model  and  the  CALSPAN  pilots'  opinion  ratings.  The  results  obtained 
from  program  OPSACT  were  compared  to  the  CALSPA”  results  in  an  attempt 
to  find  that  correlation. 

CALSPA1I  conducted  eighty-three  evaluations  using  two  flight  profiles. 
Profile  1  consisted  of  an  ILS  approach  followed  by  a  touch-and-go  landing. 
The  pilot  first  evaluated  the  overall  approach  and  landing,  and  then  he 
rated  the  approach  separately.  Profile  2  consisted  of  an  ILS  approach 
only.  Since  the  task  modeled  in  this  thesis  was  just  the  ILS  approach, 
the  overall  rating  given  for  Profile  1  flights  could  not  be  used  as  a 
basis  for  comparison  with  the  optimal  pilot  model  results. 

There  were  forty-nine  aircraft/control  system  configurations  that 
were  evaluated,  some  of  which  were  rated  using  both  flight  profiles. 

CALSPAJJ  noted  that  in  the  cases  where  a  single  configuration  was  evaluated 
using  two  profiles,  the  approach  rating  made  on  the  Profile  2  flight  was 
higher  (worse)  than  the  separate  approach  rating  made  on  the  Profile  1 
flight  (Ref.  l).  This  indicated  that  the  approach  rating  for  the  latter 
flight  profile  could  have  been  influenced  to  some  degree  by  the  landing 
instead  of  being  based  solely  un  the  Cooper-Harper  criteria.  Since  the 
model  predictions  were  designed  to  be  compared  to  the  pilot  opinion  ratings 
of  the  ILS  approach  based  only  on  Cooper-Karper  criteria,  Profile  1  approach 
ratings  ^ould  not  be  used  t:  compare  with  the  model  results.  This  left 
only  eleven  evaluations,  (eight  aircraft/control  system  configurations 
evaluated  by  two  pilots),  fcr  comparison.  The  results  of  these  evaluations 
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are  summarized  below. 


Pilot/Fit  No. 

Configuration 

Approach  Rating 

A/1883 

2-1 

3 

A/1890 

5-3 

B/1901 

5-3 

3 

A/1899 

1-6 

5 

A/1889 

2-H 

5 

B/1892 

2-1+ 

3 

A/189H 

1+-1+ 

5 

A/I885 

1-3 

6 

B/1898 

1-3 

6 

A/I887 

3-2 

6 

A/1897 

3-3 

7 

Figure  23  CALSPAN  Evaluations  Used  For  Comparison 
With  Optimal  Pilot  Model  Results 


Program  OPSACT  simulated  eight  optimal  pilot/aircraft  config¬ 
urations  performing  an  ILS  approach.  Using  CALSPAN  nominal  data,  the 
simulation  was  run  for  200  seconds  to  approximate  the  duration  of  the 
actual  approaches.  Run  1  used  arbitrarily  chosen  initial  conditions  of 
+5  ft /sec  in  forward  velocity  and  +  25  ft  in  altitude.  The  model  results 
are  listed  in  Fig.  2b.  These  results  are  then  compared  to  CALSPAN  pilot 
A's  opinion  ratings  in  Pig.  25.  The  Cooper-Harper  and  performance  index 
scales  were  aligned  by  first  performing  a  least-squares  linear  regression 
for  each  plot,  (that  is,  find  the  straight  line  that  minimizes  the  sum 
of  the  squares  of  the  deviations  of  the  actual  data  points  from  the  strai 
line  of  best  fit).  The  two  vertical  scales  were  then  aligned  so  as  to 
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superimpose  the  best  fit  straight  lines  from  both  plots.  Correlation 
between  the  pilot  ratings  and  the  model  predictions  would  be  indicated 
if  the  two  plots  were  similar  in  shape  and  slope. 


Configuration 

Performance  Index 

2-1 

7,113 

5-3 

10,110 

1-6 

8 ,66o 

2-b 

11.981 

h-h 

12,621 

1-3 

8,306 

3-2 

lk  ,266 

3-3 

13,082 

Initial  conditions: 

+5  ft/sec  (airspeed) 

+25  ft  (altitude) 

Figure  2h  Optimal  Pilot  Model  Results:  Run  ^1 


As  can  be  seen,  the  plots  seem  to  be  similar  except  for  configurations 
1-6,  3-3,  and  especially  1-3.  The  mean  difference  between  the  two  plots 
is  equivalent  to  1.1  Cooper-Harper  rating  increments.  To  determine  whether 
this  difference  between  a  human  pilot  and  the  optimal  pilo4'  model  is 
reasonable,  the  mean  difference  between  the  ratings  of  those  aircraft/ 
control  system  configurations  evaluated  by  both  CALSPAN  pilots  was  examined 
Of  the  forty-nine  configurations  flown  by  CALSPAfJ,  seventeen  were 


evaluated  by  both  pilots.  The  results  of  these  evaluations  are  depicted 
in  Figures  26a  and  26b. 


Figure  25  Comparison  of  Optimal  Pilot  Run  #1  Results 
to  CALSPAN  Pilot  A’s  Opinion  Ratings 


The  mean  difference  between  the  opinion  ratings  of  the  two  CALSPAN  pilots 
was  1.1  Cooper-Harper  rating  increments.  The  largest  difference  in  ratings 
between  the  optimal  pilot  model  and  CALSPAN  pilot  A  was  2.5  rating  incre¬ 
ments,  (configuration  1-3),  whereas  there  are  four  configurations  where 
the  two  CALSPAN  pilots  differed  by  2.5  or  more  rating  increments,  (config¬ 
urations  3-C,  3-1,  *+-3,  and  7-3).  This  indicates  that  significant  differences 
between  the  model  and  the  pilot  are  not  unreasonable.  The  fact  that  the 
mean  difference  in  ratings  for  the  model-pilot  A  comparison  and  the  pilot  B- 
pilot  A  comparison  were  the  same  could  indicate  that  the  model  is  Just  as 
accurate  in  predicting  pilot  A  ratings  as  pilot  P.  The  limited  number  of 
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Configuration 

Pilot  A  Rating 

Pilot  B  Rating 

2-1 

2 

2 

3-C 

2 

5 

1-C 

3 

3 

1-C 

1 

1 

1-1 

1 

1 

2-2 

1.5 

1 

Profile  #1  2-A 

Overall  Rating 

1 

6 

7-3 

1 

6 

3-1 

1.5 

7 

5-1 

7 

5 

5-3 

6.5 

6 

U — U 

7 

6 

3-6 

7 

6 

2-7 

7 

6 

1-3 

6 

8 

1-3 

9 

10 

1-C 

1.5 

2 

2-2 

2 

2 

2-A 

2 

3 

5-3 

3 

2 

Profile  n  1-3 

2 

5 

Approach  Rating  ^-3 

2 

6 

2-7 

1 

3 

1-1 

1.5 

3 

3-6 

5 

5 

5-3 

1 

3 

Profile  #2 

Approach  Rating  2-1 

5 

3 

1-3 

6 

6 

Figure  26a  Pilot  A  -  Pilot  B  Comparison 


Figure  26b  Pilot  A  -  Pilot  B  Graphical  Comparison 

evaluations  available,  however,  make  such  an  assumption  premature.  Also, 
the  effects  of  arbitrarily  selecting  the  initial  conditions  is  unknown. 

A  second  run  of  OPSACT  was  made,  this  time  changing  the  initial  conditions 
to  -20  ft  in  altitude.  The  results  of  Run  2  are  depicted  in  Fig.  27  and 
a  comparison  of  these  results  to  pilot  A's  ratings  is  made  in  Fig.  28. 

The  vertical  scales  were  aligned  using  the  same  method  used  for  Run  1. 

The  mean  difference  in  ratings  between  the  model  and  the  pilot  increased 
to  1.6l  Cooper-Harper  rating  increments.  Clearly,  the  use  of  mean  difference 
in  ratings  as  a  means  of  quantifying  the  degree  of  correlation  between  the 
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Configuration 

Performance  Index 

2-1 

2030 

5-3 

3483 

1-6 

2940 

2-4 

3575 

4-4 

3479 

1-3 

2761 

3-2 

3612 

3-3 

3430 

Initial  conditions: 

-20  ft  (altitude) 

Figure  27  Optimal  Pilot  Model  Results:  Run  #2 


optimal  pilot  model  and  a  human  pilot's  opinion  ratings  is  unreliable 
unless  the  effects  of  the  initial  conditions  are  known.  Unfortunately, 
the  initial  conditions  of  the  CALSPAN  evaluations  were  not  known  and 
therefore  could  not  he  used  in  the  computer  simulations.  Since  the 
initial  conditions  of  the  optimal  pilot  model  appear  to  have  a 
significant  impact  on  the  rating  predictions,  it  is  logical  to  assume 
that  the  model  initial  conditions  would  have  to  be  the  same  as  those  in 
the  actual  flight  tests  before  a  comparison  between  the  model  results 
and  the  pilot  ratings  could  be  made.  An  *r  possibility,  though,  would 
be  to  compare  the  pilot  ratings  tc  .  >me  res^t  of  the  optimal  pilot  model 
that  was  independent  of  the  initial  conditions. 

The  state  covariance  matrix  is  not  dependent  on  the  initial  conditions 
of  the  model.  It  is  dependent  only  on  the  system  conf iguraticn  and  its 
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flying  qualities.  The  RMS  error  in  forward  airspeed  ^  ,  downward  airspeed 
<7^,  pitch  rate  pitch  attitude  ,  and  altitude  <7^  were  calculated  by 
OPSACT  and  are  listed  in  Fig.  29.  Plots  of  the  errors  versus  aircraft 
configuration  are  depicted  in  Figs.  30  -  3^ .  Though  the  errors  are 
different  in  magnitude,  the  shapes  of  the  first  four  plots  are  quite 
similar.  The  plot  of  the  error  in  altitude  is  nearly  the  reverse  of 
the  other  plots.  This  may  indicate  that  the  system  that  responds  to  noise 
with  very  small  state  perturbations  will  allow  the  large  altitude  changes 
that  result  from  the  noise,  whereas  the  system  that  is  highly  responsive 
to  the  noise  will  consequentially  keep  altitude  perturbations  small. 

Figure  35  shows  a  comparison  of  the  RMS  error  in  pitch  response 
to  noise  with  CALSPAN  pilot  A’s  opinion  ratings.  The  verticle  scales  were 
aligned  in  the  same  manner  as  with  the  performance  index  comparisons.  The 
mean  difference  in  this  case  was  1  Cooper-Harper  rating  increment,  a 
difference  similar  to  that  of  the  Pilot  A  -  Pilot  B  comparison.  Also,  it 
can  he  seen  that  the  RMS  error  agrees  with  both  performance  index  comparisons 
by  predicting  a  much  lower  (better)  rating  for  Aircraft/Control  System 
Configuration  1-3  than  both  CALSPAN  pilots.  While  it  can  be  pointed  out 
that  the  two  human  pilots  disagreed  on  occasion  by  a  greater  amount  than 
the  rating  difference  for  Configuration  1-3,  and  that  even  with  this 
difference  the  RMS  error  predicts  pilot  A's  opinion  ratings  with  an  average 
error  of  1  rating  increment,  the  small  sample  size  of  this  thesis  demands 
a  low  degree  of  confidence  in  this  conclusion.  Furthermore,  a  closer  look 
at  the  comments  made  by  the  pilots  reveals  a  lesser  amount  of  correlation 
between  the  optimal  pilot  model  predictions  and  the  pilot  opinions. 

Figu.  e  36  depicts  the  same  results  as  Fig.  35  but  with  two  additions: 
CALSPAN  pilot  B's  opinion  ratings  and  a  summary  of  the  comments  made  by 
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Configuration 

<T 

V 

2-1 

11.50 

19. hh 

•  1399 

.0671 

11.29 

5-3 

11.  lU 

15. 31* 

.0397 

.OU59 

21.  U5 

1-6 

11.  lU 

15.31 

.0327 

.0550 

18.55 

2-U 
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Figure  29  HMS  Error  in  State  Response  to  System  Noise 
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Figure  30  RMS  Error  in  Forward  Velocity  Response 
to  System  iloise  versus  Aircraft/Control  System  Configuration 
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may  have  been  an  indicator  of  overcontrol  problems,  but  the  high 
RMS  error  for  Configuration  h-h  contradicts  this  assumption.  The 
problems  encountered  in  trying  to  find  correlation  between  the  optimal 
pilot  model  results  and  the  human  pilot's  opinion  ratings  (RMS  error 
compared  to  pilot  ratings)  are  again  found  in  trying  to  find  some 
connection  between  the  RMS  errors  and  the  pilots'  comments.  The  small 
sample  size  of  this  thesis  leaves  doubt  about  any  conclusion. 
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Figure  36  RMS  Error  -  Pilot  Rating  Comparison 
with  Pilots'  Comments 

A  final  attempt  to  find  some  method  of  predicting  pilot  ratings  from 
model  results  was  made  by  plotting  the  trace  of  the  matrix  resulting  fror  tie 
product  of  the  state  covariance  matrix  V  and  the  weighting  matrix  0.  No 
definite  correlation  of  any  type  was  observed. 
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V  Conclusions  and  Reeommer-Ju'  ions 

Due  to  the  small  sample  size  of  this  thesi.  ,  the  cont';  :en  e  level 
for  any  conclusion  made  would  be  extremely  low.  V.hile  some  *  rencs,  (or 
perhaps  suggestions  of  trends),  can  be  seen  in  the  results,  significant, 
contradictions  are  also  evident.  With  the  scarcity  of  data  available  to 
compare  to  the  optimal  pilot  model  predictions,  this  thesis  lias  succeeded  ; u 
primarily  pointing  to  areas  where  farther  work  vorl  is  nee  led.  !  !!,  •mi.  i  i.  ■  i.p 

have  been  learned  about  possible  methods  of  predicting  human  {  ilot  o:  inion 
ratings  based  on  optimal  pilot  model  results. 

Using  the  performance  index  as  a  measur  ■*  ot  the  a  ire  tv  ft1::  flying 
qualities  may  require  that  the  initial  conditions  of  the  tic h  1  !•  the  sar 
as  for  the  actua L  flight  tests.  This  doesn’t  s  -ein  re..  ,<>n.< l> i e  in  •  ,  I  .  i  •,(  . 
one  would  expect,  that  two  human  pilots  wool'!  give  the  cam*  rut  it.,  to 
aircraft  even  if  the  initial  conditions  of  *  he  flight  tasks  v**re  ii f feren* 
and,  second,  if  the  performance  index  were  nr.  'curate  men. 3  of  j  r*-d  1  -t  i  ug 
flying  qualities  ratings,  it  too  should  give  t  e  rune  rat  ng  as  a  human 
pilot  even  if  the  initial  conditions  were  different.  The  ‘wo  'A'  TAN  1  i  4 s 
did,  however,  vary  significantly  in  their  raMrgs  f  several  uir  tuft.  '  ••>• 
possibility  that  initial  conditions  do  have  a  irr  •  impact  -,r.  u  1  '  '  ■*  ’  s 
ratings  still  exists  utid  should  he  investigate.!  fn  -t.her. 

The  RMS  -rror  in  state  response  to  syst  ;m  1  oise  ap  eu .-*-•)  l'-;rr-  ■!  •» *  • 

with  CALSPA11  pilot  A’s  opinion  ratings  as  well  as  lid  pil  t  P’s  ■'‘tings, 
though  the  problem  of  small  sample  size  generates  doubt  ab.  ut  t>  relial:’.  :ty 
of  this  correlation.  Furthermore,  the  lack  of  cot  relation  I  '  w*-  t  the 
error  and  the  p. lots'  comments  about  the  aircraft:  '  flying  q.alities  ad  is 
to  this  doubt.  Clearly,  more  pilot  evaluat  on:  a:  e  need*.-  1  b-  for  -  anything 
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definite  can  be  concluded.  Unfortunately,  calculation  of  the  state 
covariance  matrix  (Eq.  27)  makes  assumptions  that  limit  useable  flight 
evaluations  to  those  using  long,  constant  angle  glide  slopes  without 
roundout  and  landing.  Evaluations  other  than  those  used  in  this  thesis 
cannot,  therefore,  be  used  for  comparison  with  RMS  error  in  state  response 
to  system  noise.  In  order  to  use  other  CALSPM  evaluations,  the  optimal 
pilot  model  would  have  to  be  modified  at  those  points  where  an  infinite 
glide  slope  was  assumed  and  provisions  would  have  to  be  made  to  account 
for  weighting  matrices  that  change  with  respect  to  time  and  the  forcing  of 
the  states  to  their  nominal  values  at  the  final  time.  The  model  would 
then  simulate  the  entire  ILS  approach  and  landing  and  allow  the  use  of  all 
Profile  #1  CALSPAN  evaluations  for  comparison  with  the  model  results. 

Finally,  it  should  be  noted  that  even  if  all  eighty-three  CALS  FA!.' 
evaluations  could  be  used  for  comparison  to  the  optimal  pilot  model  results, 
this  data  would  still  be  deficient  because  only  two  human  pilots  were 
used  and  their  debriefings  were  perhaps  a  bit  too  brief.  More  pilots 
flying  and  evaluating  the  same  aircraft/control  system  configurations  are 
needed  in  order  to  determine  the  mean  difference  between  human  pilot  ratings. 
This  would  indicate  how  accurate  the  model  predictions  actually  were.  The 
debriefing  of  the  pilots  would  have  to  be  much  more  detailed,  including 
questions  on  why  the  pilot  gave  an  aircraft  a  particular  rating,  other 
questions  should  be  directed  towards  the  aircraft  response  that  results 
from  the  RMS  error  in  state  response  to  system  noise.  Clearly,  the  use 
of  experimental  data  not  gathered  in  conjunction  with  the  computer 
simulations  adds  many  sources  of  error  that  can  completely  distort  the 
results  of  a  study.  Still,  while  this  thesis  did  not  show  that  the  optimal 
pilot  model  could  predict  pilot  ratings,  it  did  not  disprove  this  theory  either. 
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Appendix  A 


First  Necessary  Conditions  for  Optimality 


The  purpose  of  an  optimal  control  system  is  to  find  the  input  u*(t) 
that  minimizes  the  performance  index 


J(u)  =  f1[x(tf)]  + 


f2Wt),u(t),t]dt 


(A-l) 


subject  to  the  constraining  equation 


x(t)  =  f^txjt) ,u(t) ,t] 


x(t  ) 
—  o 


X 

~o 


(A-2) 


One  method  for  developing  the  First  Necessary  Conditions  for  Optimality 
is  described  by  Silverthorn  (Ref.  9).  The  state  variables  x(t)  and  the 
control  inputs  u(t)  are  treated  as  being  independent  by  introducing  a 
vector  of  Lagrange  multipliers  X(t)  to  augment  the  performance  index 
with  the  constraining  equation 


J(u)  =  f1[x(tf)]  + 


,u(t )  ,t  ] 


+  XT(t) [ f  [x(t ) ,u(t ) ,t ]  -  x(t)] 


(A-3) 


The  Hamiltonian  is  defined  by 
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H[x(t),u(t),t]  =  f2[x(t),u(t),t] 

+  XT(  t )[  f^[x(t)  ,u(t )  ,t]  ]  ( A— 1+ ) 


Equation  A-3  becomes 

J(u)  =  f1[x('tf)] 

■tf 

+  J"  H[x(t )  ,u(t )  ,A(t )  ,t ]  -  AT(t)x(t)|  dt  C A— 5 ) 

t 

o 

Recalling  that  u*(t)  is  defined  as  the  control  input  that  minimizes 
J(u),  then  x*(t)  and  £*{t)  represent  the  corresponding  optimal  trajectories. 
Non-optimum  values  of  these  variables  separated  from  the  optimum  by  a  small 
variation  are  represented  by 


u(t)  =  u*(t) 

+  €l(t) 

(A-6) 

x(t)  =  x*(t) 

+  6  £(t) 

(A-7) 

A(t)  =  A*(t) 

+  €  »(t) 

(A-8) 

Substituting  the  relations  in  Eqs.  A-6,  A-7,  and  A-8  into  Eq.  A-5  and 
expanding  in  a  Taylor  series  yields  Eq.  A-9  on  the  following  page.  The 
First  Necessary  Conditions  for  Optimality  are  based  on  the  following 
argument.  By  choosing  €  arbitrarily  small,  the  C  term  in  Eq.  A-9,  unless 
it  is  zero,  will  overwhelm  the  higher  order  terms.  Since  £  can  be  chosen 
either  positive  or  negative,  J(u)  evaluated  at  some  non-optimum  u(t)  will 
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be  greater  than  cr  equal  to  J(u*)  only  if  the  coefficient  of  the  tern  5s 
zero.  Thus,  for  u*(t)  to  e  an  optimal  control,  the  following  must  he  true 


-  X*T(t)i(t)  -  2tr(t)i*(t) 


0 


( A-10) 


The  values  f_(t)  and  £(t)  cannot  be  chosen  independently,  so  by  integrating 

* 

the  ?_(t)  term  by  parts  and  rearranging,  Eq.  A-10  becomes 


(&) 


'>*4 


-  x»T(tf)U(tf)  + 


+ 


)j>£(t)dt  + 


( A-ll ) 


Since  the  selection  of  £(t)  is  arbitrary,  one  choice  is  to  select  £(t^)  =  0. 
Each  of  the  remaining  integrals  in  Eq.  A-ll  must  be  equal  to  zero. 

Applying  the  Fundamental  Lemma  of  the  Calculus  of  Variations  to  the 


integrals  yields 
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(A-12) 


(A-13) 


(A-lU) 


Another  choice  for  i_(t)  is  ^(t^)  #  0.  Thus,  the  first  term  in  Eq.  A-ll 
becomes 


(A-15) 


Also,  recall  that 


x(t  )  =  x  (A-l6) 

—  o  —o 

Equations  A-12,  A-13,  A-lU,  A-15,  and  A-l 6  are  the  First  Necessary 
Conditions  for  Optimality.  The  general  optimal  control  problem  could  now 
be  solved,  except  that  the  boundary  conditions  represented  by  Eqs.  A-15 
and  A-l6  are  "split",  (one  boundary  condition  is  given  at  the  initial 
time  and  the  other  at  the  final  time).  The  solution  must  be  guessed  at, 
compared  to  the  First  Necessary  Conditions  for  Optimality,  and  then 
iteratively  reguessed. 

The  exception  to  the  iterative  solution  process  occurs  when  the 
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t 


system  is  linear  and  the  performance  index  is  a  quadratic  form.  The 
general  form  of  the  linear  system,  quadratic  performance  index  is 


J(u)  =  \  x" (t  )Hx(t  ) 
f  -  f 


+ 


'2 


[x1(t)Qx(t)  +  u^(t )Ru(t ) ]dt 


(A-17 ) 


The  Hamiltonian  is 

H[x(t) ,u(t) ,A(t) ,t]  =  \  xT(t)Qx(t) 

+  %  uT(t)Ru(t)  +  AT(t)[Ax(t)  +  Bu  ( t ) ]  (A-l8) 

where  the  constraining  equation  is 

x(t)  =  Ax(t)  +  Bu(t)  ;  x(t  )  =  x  (A-1Q) 

—  —  —  —o—o 

Applying  the  First  Necessary  Conditions  for  Optimality  yields 

A*(t)  =  -  ^  =  -  Qx*(t)  -  ATA*(t)  (A-20) 

*2.1# 

0  =  Ru* ( t )  +  BTA*(t)  ( A-21 ) 

x(t)  =  Ax(t )  +  Bu ( t )  (A-22) 

A*(tf)  -  Hx*(tf) 
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x(t  )  =  x 

—  O  — D 


(A-2U) 


Eq.  A-21  becomes 


u*(t)  =  -R_1BTA*(t)  (A-25) 

Combining  Eqs.  A-20,  A-22,  and  A-25  into  matrix  form  yields 


x(t) 

r  •  -1  t  i 

A  1  -BR  B 

1 

x(t) 

i(t) 

# 

1  rn 

-Q  1  -A 
-  1  « 

x{t)  * 

The  expression  for  A(t)  is  assumed  to  be  a  linear  combination  of  x(t); 
thus , 


A(t)  =  K(t)x(t) 


(A-2T) 


Eq.  A-26  becomes 


x(t) 

1  -1  T 

A  1  -BR  B 

1 

x(t) 

K(t)x(t)  +  K(t)x(t) 

V  f 

_-Q  .  -Ax 

_K(t)x(t) 

(A-28) 


Reducing  Eq.  A-28  to  a  single  equation  yields 


-K(t)  =  K(t)A  +  ATK(t)  +  Q  -  K(t)BR~1BTK(t) 


(A-29) 
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£qs.  A-27  and  A-23  imply  that 


K(tf)  =  H 


(A-30) 


and  combining  Eqs.  A-25  and  A-27  yield 


u*(t)  =  -R_1B'rK(t)x(t) 


(A-31) 


If  the  final  time  is  far  off  into  the  future  relative  to  the  time  constants 
of  the  system,  K(t)  can  be  considered  a  constant.  Thus,  the  weighting  cn 
the  final  states  is  irrelevant  and  Eq.  A-29  becomes 

0  =  KA  +  ATK  +  Q  -  KbF-1BTK  (A-32) 


and  Eq.  A-31  becomes 


u*(t)  =  -R_1BTKx(t) 


( A— 33 ) 
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Appendix  B 


Mathematical  Development  of  the  Estimator  and  the  Predictor 


As  mentioned  in  Chapter  II,  the  optimal  pilot  model  evolves 
from  the  optimal  control  system  through  the  addition  of  an  estimator 
and  a  predictor.  The  mathematical  development  of  these  two  components 
is  now  discussed. 

Estimator 

Since  the  system  states  cannot  he  directly  observed  by  the  pilot, 
they  must  be  reconstructed  in  his  mind  based  on  some  known  information. 
His  estimation  of  what  the  states  are  is  modeled  using  a  full-state 
observer.  The  observer  reconstructs  the  system  states  x{t) ,  which 
are  defined  by 


x(t)  =  A  x(t)  +  B  u ( t )  ;  x(t  )  =  x  (B-l) 

—  1—  1—  —  o  — o 

and  produces  an  estimate  of  the  states  x(t)  according  to 

x(t)  =  A  x(t)  +  B  u(t)  ;  x(t  )  =  x  (B-2) 

—  1—  1—  —  o  — o 

Using  the  same  input  u(t)  as  the  system  and  setting  x  =  x  ,  the 

—  — o  — o 

observer  should  provide  an  exact  estimate  x(t)  =  2l(t).  However,  due  to 
noise,  modeling  inaccuracies,  and  initial  conditions  an  exact  estimate 
is  not  possible  (Ref.  13).  To  improve  the  estimate,  more  information 
about  the  system  is  available  to  the  pilot;  namely,  the  displayed  variables. 


Recall  that 
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£(t)  =  Cx(t ) 


The  estimate  of  ^(t)  is 


A,  .  A 

2.(t)  =  Cx(tl 


(B-3) 


(B-U) 


The  error  in  estimating  the  displayed  variables  is 

^(t)  =  i(t)  -  £(t) 


(B-5) 


To  account  for  modeling  errors  and  noise,  a  correction  term  is  added 
to  Eq.  B-2 


x(t)  =  Ajc(t)  +  B^u(t )  +  S[^e(t)]  (B-6) 

where  S  is  the  estimator  gain  matrix  and  will  be  discussed  in  more  detail 
later.  By  defining  the  error  in  estimating  the  states  as 

x^t)  =  x(t)  -  x(t)  (B-7) 

it  can  be  seen  that 

x^Ct)  =  x(t)  -  x(t)  =  A^Ct)  -  S[}^(t)] 

=  (A, -SC )x  (t)  (B-8) 
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or  equivalently 


(A-SC)(t-t  ) 

x  (t)  =  e  x  (t  )  (B-9) 

- e  —eo 

The  estimate  error  will  decay  to  zero  only  if  S  is  chosen  such  that  the 
system  described  by  Eq.  B-8  is  stable  (Refs.  13,11*).  Assuming  complete 
observ-bility  of  the  system,  the  S  matrix  can  be  chosen  such  that  the 
observer  poles  defined  by  the  eigenvalues  of  the  matrix  [A^-SC]  are  far 
enough  into  the  left-hand  complex  plane  to  insure  that  the  decay  of  the 
error  is  fast  enough,  yet  not  so  far  to  the  left  as  to  make  the  S  matrix 
large  and  the  observer  very  sensitive  to  observation  noise  (Ref.  lU). 

Before  determining  the  gain  matrix  S,  it  should  be  recalled 
from  Chapter  II  that  the  system  is  actually  a  stochastic  process.  Random 
variables  are  introduced  into  the  system;  therefore,  Eqs.  B-l  and  B-3 
are  rewritten  as 


c(t)  =  A^x ( t )  +  B^u ( t )  +  D2vg  (t) 


(B-10) 


jr(t)  =  Cx(t)  +  ]Mt) 


(B-ll) 


where  w  (t)  and  D  are  as  described  by  Eq.  2k,  and  v  (t)  represents 
observation  noise.  Both  processes  are  assumed  to  be  zero-mean,  white 
Gaussian  noise  of  intensities  and  V^,  respectively.  The  optimal 
observer,  as  defined  by  Eq.  B-6,  attempts  to  minimize  the  mean  square 
error 


E(x  ^(t)Gx^(t) ^  =  x^(t)  Gx^(t)  +  tr[PG] 


(B-12) 
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where  G  is  a  symmetric,  positive  definite  weighting  matrix,  x^Ct)  is  the 
mean  of  x  (t),  and  P  is  the  estimator  error  covariance  matrix  (Ref.  li*). 


P(t) 


=  E{ 


[x  (t)  - 
— e 


X  (t)][x  (t)-x  (t)]' 
— e  -e  — e 


(B-13) 


As  mentioned  before,  apart  from  the  effects  of  noise,  model  inaccuracies, 
etc.,  the  estimator  error  can  he  minimized  by  the  proper  choice  of  the 
initial  conditions.  The  first  term  of  Eq.  B-12  is  minimized  when 
x^  =  x^  The  stochastic  inputs  in  Eqs.  B-10  and  B-ll  change  the  estimator 
error  described  by  Eq.  B-8  to 


x  (t)  =  (A, -SC)x  (t)  +  D0w  (t)  -  Sv  (t) 
~e  1  -e  2  g1  -y 


=  (A-SC)Z(t)  +  F 

X  “c! 


v  (t ) 

y 


(B-lU) 


where  F  =  jD^  J  -S 

In  this  study,  the  estimator  is  assumed  to  have  been  operating  for 
some  time  beyond  the  start  of  the  flight.  Therefore,  it  can  be  assumed 
that  the  estimator  error  covariance  matrix  P  has  reached  a  steady  state 
value,  (much  the  same  way  the  solution  to  the  optimal  control  Riccati 
equation  K  can  be  considered  to  be  constant  if  the  final  time  is  far 
into  the  future).  Thus,  it  can  be  shown  (Ref.  ll*)  that  the  s+eady  state 
estimator  error  covariance  matrix  is 


8U 
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“  (A  -SC)t  (A  -SC)Tt 

P  =  I  e  1  FW2?  e  1  dt,  w  = 


ot 

■J 


w  Jo 

-*•  i 

—  T. — 

0  •  V 

:  y 


(B-15) 


The  covariance  matrix  P  satisfies  the  equation 


0  =  (An-SC)P  +  P(A  -SC)  +  FW  v 
1  1  2 


=  (A,  -SC)P  +  P{A,  -SC)T  +  W,  +  SV  sT 
1  1  1  y 


(B-l6) 


Furthermore,  P  is  minimized  if  (Bef.  ll) 


T  -i 
S  =  PC  V 

y 


(B-1T) 


Substituting  Eq.  B-17  into  Eq.  B-l6  yields 


0  =  j[P  +  PA^  +  W1  -  PCTV"1CP 


(B-l8) 


The  value  of  P,  when  applied  to  Eq.  B-17,  yields  an  estimator  gain  matrix 
S  that  minimizes  the  second  term  in  Eq.  B-12.  Figure  37  shows  how  the 
estimator  fits  into  the  optimal  pilot  model. 

Predictor 

When  a  delay  T  exists  in  the  system,  the  estimator  produces  the 
delayed  estimate  of  the  states  x(t-T),  instead  of  |(t).  Eq.  B-6  becomes 


x(t-T)  =  A^t-T)  +  B^u( t-T )  +  Sl^(t-q')  -  Cx(t-T) ] 


(B-19) 
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y(t)  =  Cx(t) 


Optimal  Pilot  Model 
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The  predictor  models  the  pilot's  reconstruction  of  the  system  states  x(t) 
based  on  the  estimate  of  the  delayed  states  x(t-T).  The  predictor  is 
developed  by  Kleinrcan  (Ref.  ll)  by  first  defining  x*(t)  as  the  part  of 
x_(t)  due  to  the  deterministic  input  u*(t),  where  x_*(t)  satisfies  the 
relation 


x*(t)  =  A  x*(t)  +  B  u*(t) 


(B-20) 


The  predictor  error  is  then  defined  as 


r(t)  =  x(t)  -  x(t) 


(B-21) 


which  is  a  random  process  and  evolves  according  to 


r(t)  =  A^r(t)  +  S[jr(t)-Cx(t)] 


(B-22) 


Eq.  21  can  be  written  as 


x(t)  =  x*(t)  +  r(t) 


=  x*(t)  +  E  |r(t)|  x(tf) ,  O’ it  | 

=  x*(t)  +  E  jr(t)  x(  O'  ),  Ci  t-lj> 


(B-23) 


Since  r_(t)  is  independent  of  u*(t),  Eq.  B~23  becomes 


=  x*(t)  +  E  jr(t)  r(  ^ 


=  x*(t)  +  e  r(t-T) 


(B~ph) 
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The  reconstructed  state  estimate  is  determined  fcy  Eq.  B-2h ,  Figure  37 
shows  how  the  predictor  is  added  to  the  optimal  pilot  model. 
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Appendix  C 


-  Single  Axis  Control  Task  (OPSACT)  Profirar  Listing 


Optimal  Pilot 
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1 

2 
3 

s 


5 

b 

7 

a 

9 

10 

U 

12 

13 

1* 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2*. 

25 

26 

27 

28 

29 

30 

31 

32 

33 
35 

35 

36 

37 

38 

39 
80 
81 
82 
83 
88 

85 

86 

87 

88 
89 

50 

51 

52 

53 
58 

55 

56 

57 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 


c 

c 

c 

c 

c 


PROGRA*  OPS  AC  T  < 1 NPUT-/80. OUTPUT ,PL OT , TAPES- INPUT ,T APEb-OUTPuT) 

PROGRA-  OPSACT  (OPTI-AL  PILOT  -  SINGLE  AXIS  C  ON  7  POL  TASK), 

AT  THIS  STACE  op  its  oe velop^ent ,  IS  DESIGNED  TO  "OOEl  THE 
OPTI-Al  PILOT  PEPEOR“ING  AN  INFINITE  T  1 78  £  LINEAR  REGULATOR  TASK. 
THE  input  values  are: 

III  SYST£“  COEFFICIENT  8ATRIX  ( AO > 

(21  CONTROL  INPUT  8  A  T  R It  (80) 

(31  OUTPUT  8f  A  SURE"  ENT  8A  TR  I  x  |C> 

(81  STATE  WEIGHTING  -ATRfX  (00) 

(51  CONTROL  WEIGHTING  8ATRIX  (RO) 

(6)  OBSERVATION  nOISE  COVARIANCE  "ATPJX.(V) 

(71  EXTERNAL  DISTURBANCE  COVARIANCE  8ATRIX  (WO) 

(8  1  INITIAL  T I 8£  ( T I  ) 

(9|  final  T  I  8£  (TF) 

(101  INTEGRATION  T I 8£  INCRE-ENT  ITINCR) 

(111  INITIAL  STATE  VALUES  (XI) 

(121  INITIAL  CONTROL  INPUT  VALUES  IUI) 

(13)  INSTRU-ENT  THRESHOLDS 

THE  SrSTE*  IS  PRESENTLY  L 1 8 1  TED  TO  A  TWELFTH  ORDER  PROBLE  8 . 
REAL  A(I2.12I.B(12,121,C(12.12).0I12.12).V(12.12),R(12,12), 

W(12.l2).P(12.12).K(12,12),sri2tl2),l(12.12),PHU12,l21. 

0REV(12,12I,GA88A(12,12),PSI(12.12).RESULTI12,12), 
RtNV(12,12).8TRANS(12.12).RESULT2(12.12) .COUTI 12 .1 2) , 

CTR ANSI  12.121 .VINV( 12.1 2) ,ATR ANSI  12,12 )  .TAUNI 12,12) , 
U2(12.12),X2(12,12),2I12,12).21I12.12).RESULT3(12,12), 
XOELAVI 12.12) *  UOEL AT I 12 . 1 2  > . TOE L A T I  1 2 , 1 2  > , 

XNEW( 12) .UOLOI  12 ).OXHATNI 12 ) .DXHATOI 12 )  .7(12 ) ,  7 
TEST! 12  I .OEE ( 12 ) .EE( 12 ) .XHATI 12 ) .OXHATNI 12) ,UQPT( 12 ) , 

OUTPUT! 12) ,72 ( 12), 23(12), 78112), THRSHD 112), 

X0L0( 12  I ,0XHAT0( 12 ) , X  I  ( 12 ) ,UI ( 12) ,TP( 300) ,TA( 1200 ) , 

00XHAT ( 12  I ,OUOPT( 12 ) .YNEWI12), 

TOLO.riNCR.  TNEW.TF,  TOL,  TAU,  T  I  ,  0  T  ,  U3  ,  X  3  ,  P  F  8,  PF  «  T  ,  SU8 1 ,  SU8  2 

INTEGER  «OOE . N, «,rN,IER,J, «R, IRVING, CHANGE, I. II. INDEX, NLEN, 

:  NOUT, COUNT, 80RE,J2,LA,INC.!D.IE,IG.NP0INT5.L18IT,RUN, 

:  EIG, GUESS. “8, II  I 

character  handle*? 

CO 8  WON/ N A  I N 1 / NO I  8 , N 0 1  8 1 , 0U8 (  188) 

C0  8  80N/*AIN2/ST0RE  I  188) 

C088QN/8AIN3/EXTRA(  188) 

C0880N/ INOU/K IN.KOUT.KPU 

C08«0N/C,RAFS/T  ITLE  (  3  I.  Yl  TITLE!  2)  .Y2TITLE  (2)  ,Y3TI  TLE  t  21  «Y8T  ITLE  12) 
CO'YON/BOWR/EIG 


CHANGE-2 

*****  800E-1  IE  ON  INTEPC08;  OTHERWISE,  "ODE-2. 


PRINT* (/T2.2A, II, A/I* 
REAO*.  "OOE 


*  WE  L  CO  8E  TO  OPSACT. 
•NODE  BY  ENTERING  A 


INDICATE  INTERC08  *, 
■Mi*** 
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38  PR  tNT'( /T2,  2A  .// lbf  T2  ,A/)>‘  . 

39  :  mm  order  that  the  motor  noise  be  properly  augmented  into  •» 

60  :  ‘THE  •*«*•  »ATRIX«  THE  "A"  MATRIX  *U$T  BE  ARRANGED  AS  F  QLLO -S : ’  » 

61  :  •  •  **  .  .  »•♦*. 

62  *  •  **• 

63  S  *  •  BASIC  .  .  •*« 

66  :  •  ♦  SYSTEN  .  .  »*. 

63  *  •  *  ’  • 

66  ■  '  * .  ♦  *’* 

67  *  *  •  *  • 

68  :  •  ♦  AUGMENTING  .  . 

69  •  ♦  CONTROL  INPUTS  .  ♦*. 

70  :  •  ♦  .  •  ’ . 

71  *•  *  * . 

72  :  •  ♦  «M 

73  :  •  *  auG“ENTING  *M 

76  •  ♦  CONTROL  SYSTEM  •’  • 

73  t  *  * 

76  ;  •  ♦«* 

77  C 

78  C  *****  START  ***** 

79  C 

80  l  JA80RT-0 

81  f  f  TMf70E.fO.il  THEN  * 

82  PR INT • I/T2 .A) • . ‘ENTER  R UN/CONE  I  OUR  AT  I  ON  NUN  BE  R  TINT  EGER ) * 

83  PRINT*  1/T2.2A/T2.2M’  .‘ENTER  SYSTEM  ORDER  IN),  CONTROL  *. 

86  :  ‘ORDER  IN)  *, ‘CONTROL  SYSTEN  ORDER  I  MM  1  ,  *, 

83  :  ’and  OUTPUT  MEASUREMENT  ORDER  IYN).‘ 

86  PR!NTMT2.A/*.tT6,A,ll,AMl‘,‘ENTER  KTRAJ  AS  FOLLOWS:*. 

87  :  ‘KTRAJ  -  ‘,0.‘  If  NO  PLOT  IS  DESIRED.*. 

aa  :  ‘ktraj  ■  ‘,i.‘  if  only  a  plot  is  desireo,*, 

89  :  ‘ktraj  •  ‘,2,*  If  only  A  LISTING  IS  OESIREO,*. 

90  :  ‘KTRAJ  «  ’  ,3 .  ‘  IF  BOTH  PLOT  AND  LISTING  ARE  OESIREO.' 

91  PfllNT‘(/T2,NA,Il.A)‘.‘F0R  LISTING  OF  EIGENVALUES  AND 

92  :  ‘EIGENVECTORS  OF  THE  CONTROLLABILITY/’, 

93  ;  ‘OBSERVABILITY  GRannIan  FRO*  m*iC.  *, 

96  *.  ‘SET  "EIG"  - 

03  END  IF 

96  READ*. RUN, N.M.MM.YN, KTRAJ, EIG 

97  PRINTM//T2«A,12//I'»* RUN/CONf (GURATION  NUMBER  IS:  RUN 

9B  HO  I'- l 2 

99  NOINl-NOI M*l 

100  KOUT-6 

101  KPU-A 

102  KIN-3 

103  JJ-0 

10*.  2  IF  (CHANGE. EO.  1)  THEN 

103  if  «NonE.Eo.n  then 

106  IF  IJJ.EO.O)  THEN 

107  PRINT*«‘tF  YOU  WISH  to  CHANGE  ANY  of  the  FOLLOWING  *. 

108  :  ‘INPUTS.  Enter  the  nu“8ER  in  PARENTHESES  beside  *. 

109  :  ‘THE  INPUT.  ENTER  A  ZERO  If  NO  "ORE  CHANGES  ARE  *, 

HO  :  ‘OESIREO.’ 

111  PRINT‘(/12(T2.A,I2,A>>‘, 

112  :  ‘ABATRIX . I ‘ .1. ‘ )‘ . 

113  :  ‘8  MATRIX . ( • ,2, • I*  , 

116  *.  ‘C  NATRIX . I’  .3.’  )*  . 


GAE/AA/SCD-3 


* 


MS 

116 

117 

119 

119 

120 
121 
122 
123 
129 
12  9 
126 

127 

128 

129 

130 

131 

132 

133 

138 

135 

136 

137 
135 

139 
180 
181 
182 
193 
198 

185 

186 

187 

188 
189 

150 

151 

152 

153 
158 

155 

156 

157 

158 

159 

160 
161 
162 
163 
168 

165 

166 

167 

168 
169 
1  70 
171 


•  0  MAT* 1 X. . 

..  !'.8.'l'. 

•R  MATRIX.. 

. .1 '  ,5.'  I*  . 

' V  MATRIX.. 

..( ' .6.' 1*  . 

•W  MATRIX.. 

.  .(  '  ,  7,  •  1  •  , 

' INITIAL  ti 

E . ( * .8, ' 1  '  , 

•final  TIME 

.. (  '  ,9. '  1  »  , 

• TINCR . 

.  .  (  '  ,10,  '  1  ' 

*X( [NIT [ALI 

.  .  (  ' . 11  ,  '  )  ' 

•U( INITIAL) 

.  .  (  '  .12 , '  I  ' 

• THRE  SHQLD. 

..1 '.13,')' 

'LAST  ENTRY 

..  1  '  .0.  *  )  ' 

CO  TO  b 
ENOIF 

PR INT • ( /T2.A/  )*  .  * ENTER  NUMBER  FOR  NEXT  DESIRED  CHANCE.* 
ENOIF 

6  REAO*.  IRVING 

J  J*1 

IF  ( IRVING. LT. 1. OR. IRVING. GT. 13)  THEN 
CO  TO  135 
ENOIF 

CO  TO  I  10.20.  30.80. 50,60, 70.80,90, 100, 1 10,120.  1301 ,  IRVINC 
ENOIF 

c 

10  IF  IND0E.E0.1I  THEN  * 

PRINT*.  ‘ENTER  SYSTEM  COEFFICIENT  MATRIX  CAI.  ENTER  ROW,  *. 
:  'COLUMN,  AND  VALUE;  EACH  SEPARATEO  8T  A  COMMA. * 

ENOIF 

HANOLE*'  "A"  * 

CALL  l NP MT X( HANOLE, A, N.N. A, CHANGE. MODE ) 

IF  (CHANCE.Ea.il  THEN 
CO  TO  2 
ENOIF 

20  IF  (MDOE.EO.ll  THEN 

PRINT*,  'ENTER  CONTROL  INPUT  MATRIX  181* 

ENOIF 

HANDLE  * '  "8“  * 

CALL  INPMTXI HANOLE . 8. N.N.B, CHANCE. NODE  1 
IF  (CHANGE. EO. 1 1  THEN 
GO  TO  2 
ENOIF 

30  IF  (*nOE.E0.11  THEN 

PRINT*.  'ENTER  OUTPUT  MEASUREMENT  MATRIX  (Cl* 

ENOIF 

handle*'  "C"  * 

CALL  I NP*TX(H ANGLE ,C.VN,M, C. CHANGE ,HOOE I 
IF  (CHANGE. EO. 1)  THEN 
GO  TO  2 
ENOIF 

80  IF  (MOnF.EO.il  THEN 

PRINT*. 'ENTER  STATE  WEIGHTING  MATRIX  (01* 

ENOIF 

HANDLE*'  "O"  * 

call  i np*tx( handle . o.n.n.q, change, mooei 

IF  (CHANGE. EO. II  THEN 
GO  TO  2 
ENOIF 

50  IF  (MOOE.EO.il  THEN 
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17? 

171 

174 

175 

176 

177 
1  78 

179 

180 
181 
18? 
181 

184 

185 

186 

187 

188 

189 

190 

191 
19? 
191 

194 

195 

196 

197 

198 

199 

200 
201 
202 
201 

204 

205 

206 

207 

208 

209 

210 
211 
212 
211 

214 

215 

216 

217 

218 

219 

220 
221 
222 
221 

224 

225 

226 

227 

228 


PRINT* , "ENTER  CONTROL  WEIGHTING  MATRIX  t «  I  • 

E  NO  I  F 

MANGLE ■ •  "R"  • 

CALL  INPwtx(HANOLE , R , 1, W,R, change  »800E 1 
IF  (CHANGE .EO. 11  THEN 
GO  TO  2 
ENOIF 

60  IF  MOOE.EO.il  THEN 

PRINT*, "ENTER  OBSERVATION  NOISE  COVARIANCE  -ATRIX  (VI" 

ENOIF 

HANOLE-"  " V"  • 

CALL  INPWTX! HANDLE, V.YN.VN, V, CHANGE  MOO E) 

IF  (CHANGE .EO. 11  THEN 
GO  TO  2 
ENOIF 

70  IF  MOOE.EO.l)  THEN 

PRINT*, "ENTER  EXTERNAL  DISTURBANCE  COVARIANCE  "ATRIX  |W»" 
ENOIF 

HAMOLE*"  "W"  » 

CALL  INP*"TX(  HANDLE  ,W,N,N,W,  CHANGE  MODE  1 
IF  (CHANGE.EO.il  THEN 
GO  TO  2 
ENOIF 

80  IF  (NOOE.EO.il  THEN 

PRINT*,  "ENTER  INITIAL  TIME  (TI)  AS  A  REAL  NU-BER" 

ENOIF 
READ*,  Tt 

PRINT"! /T2.A.F6.2/I",  "THE  INITIAL  TIME  IS:  ",  TI 
IF  (CHANGE. EO. II  THEN 

GO  TO  2  ' 

ENOIF 

90  IF  (WOOE.EO.il  THEN 

PRINT*.  "ENTER  FINAL  TME  (TFI  AS  A  REAL  NUMBER  " 

ENOIF 
R  E  AO* ,  TF 

PRINT*(/T2,A,F6.2/|","THE  FINAL  TME  IS:  *  ,  TF 
IF  (CHANGE .EO. II  THEN 
GO  TO  2 
ENOIF 

100  IF  (WOOE.EO.il  THEN 

PRINT*, "ENTER  THE  INTEGRATION  TIWE  INCREMENT  (TINCRI". 

:•  AS  A  REAL  NUMBER" 

ENOIF 

R  E  AO*  »  TI NCR 

PR INT"( 7T2, A,F6. 3/1  ", "THE  TINE  INCREMENT  IS:  * ,T I  NCR 
IF  (CHANGF.E0.1I  THEN 
GO  TO  2 
ENOIF 

110  IF  (wn0E.E0.il  THEN 

PRINT*, "FNTER  INITIAL  STATE  VALUES" 

PRINT • ( T2  ,  A  ,  12, A, A/ 1 " , " ENTER  AS  FOLLOWS:  ROW.  ",l,"  VALUE;", 
5  "  EACH  SEPARATED  BY  A  COMMA." 

ENOIF 

HANDLE-"  XI  • 

CALL  INPwTX!HANDL£,Xl,N,l,X! .CHANGE . MODE  I 
IF  (CHANGE. EO. 11  THEN 
CO  TO  2 
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229 

ENOIF 

23  0 

120 

I E  COOE.EO.il  THEN 

211 

POINT*,  'ENTER  INITIAL  CONTROL  INPUT  VALUES' 

232 

POINT* <T2>a,I2,a,a/>' , 'ENTER  AS  FOLLOWS:  ROW.’.l,' 

233 

:  '  EACH  SEPARATED  BY  COM M A  S • ' 

23* 

E  NO  I  E 

2  31 

HANCLE-*  Ut  ' 

236 

call  !n»-t* ihanole.ui .  m.i  ,ui  .change  code i 

237 

IE  (C“ANCE.£0. 1 1  THEN 

238 

CO  TO  2 

239 

EnOIE 

2*0 

130 

if  <’•onE.eo.li  then 

2*1 

PRINT*.  'ENTER  INSTRUMENT  THRESHOLDS.' 

2*2 

PRINT' <T2, A, 12. A, A/)'. 'ENTER  AS  FOLLOWS:  ROW.'.l,' 

2*3 

:  '  EACH  SEPARATED  BY  CO"«AS.' 

2** 

ENOIE 

2*1 

HANOLE«'THRSWLO* 

Zst 

Call  inp* tv ih anole , thrshd .yn.i . THRSHD .chance , "ODE  1 

2*7 

IF  (CHANGE. E0.1J  THEN 

2*8 

CO  TO  2 

2*9 

ENOIF 

210 

C 

211 

c 

. 

212 

c 

♦•••♦  CALCULATION  OF  OPTIMAL  FEEDBACK  CAIN  MATRIX  IH 

213 

c 

21* 

135 

CALL  TRANSP<9.BTRANS.N,M> 

211 

CALL  I NPNTXI HANDLE .R.m.M.RESULT .3.M00E I 

216 

CALL  CMINVI  m.m, RESULT  ,RINV,.MR,0! 

217 

ie  <«r.ne.w|  then 

213 

PRINT* </T2,A/»',  *R  NOT  FULL  RANK.*  v 

219 

IF-COOE.f 0.H  THEN 

260 

CO  TO  l 

261 

EnOIF 

262 

JA«0RT»1 

26  3 

CO  TO  210 

26* 

ENOIF 

261 

CALL  *“UL(R tNV,3TRANS , M, 8, N, RE  SUL T 1 

266 

CALL  wi/LIB  .RESULT.  N, m,N.RESULT2> 

267 

IFR-0 

268 

CALL  «»IC<N,A,RESULT2,0,K,2,IER» 

269 

IF  ( IER.EO.il  THEN 

270 

IF  (MOOE.EO.il  THEN 

271 

cn  TO  1 

272 

ENOIF 

273 

JA80PT*1 

27* 

GO  TO  210 

271 

ENOIF 

276 

PRINT'!  7T2, A/ 1 • ,  'THE  K  MATRIX  Is:' 

277 

DO  200  I-l.N 

278 

PRINT' ( JO! 3X.E 10. *1 1 (M I • Jl . J-l.NI 

279 

200 

CONTINUE 

280 

CALL  ""UL (RESULT, K, m.n.N, LI 

281 

PRINT'(10(3X,E10.*1)'.ILII.J1,J»1,N) 

282 

201 

CONTINUE 

283 

C 

23* 

C 

*****  DETERMINING  NEUROMUSCULAR  LAGS.  TAUINl  ***** 

281 

C 

value;' . 


value:'  . 
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’86 

11*1 

4.87 

00  203  4- l I , IN-H*) 

288 

00  202  J*  I  8-11 , «  N-«N) 

289 

T  AUNI I , Jl-l/L 1 1 ! , J) 

290 

202 

CONTINUE 

291 

I  I  -  1  I  ♦  l 

29  2 

203 

CONTINUE 

293 

C 

294 

C 

*****  FORCING  A1  AND  81 

2  96 

C 

296 

it*i 

297 

00  206  1« | N-N«)-l M-l ) , ( N-X* ) 

298 

00  205  J*l N-N4I-CN-1) , IN-NN) 

299 

Alt, 31--1/TAUN1 1 , J  » 

300 

91  I,  m-i/TAUNI  1,JI 

301 

11*11*1 

302 

205 

CONTINUE 

303 

II- 1 

304 

206 

CONTINUE 

308 

C 

306 

PRINT*I//T2,A) *,*A  (ONE)* 

307 

00  207  1*1, N 

308 

PRINT* ( 101 3X,E 10.41) •, 1 Al 1 , J1 , J-l.NI 

309 

207 

CONTINUE 

310 

PRINT* I//T2.A) *,*8  1  ONE ) * 

311 

00  209  1*1, N 

312 

PRINT*! 101 3X.E10.4)) *,191 1,31 , J*1 ,N) 

313 

208 

CONTINUE 

314 

C 

315 

c 

*****  FORNING  10  «♦♦**  v 

316 

c  _ 

317 

00  210  !»l,H 

318 

00  209  J»IN-NN)-IN— l)»l N— Ml | 

319 

HI  ,J  1-0. 

320 

209 

CONTINUE 

321 

210 

CONTINUE 

32  2 

PR  INT* 1 7T2,A1 ’ ,*  THE  L  IZEROl  MATRIX* 

323 

00  211  1*1,8 

324 

PRINT’ 1 101 3X.E10. 4 > ) » , <L I  I . J) . J-l.N) 

325 

211 

CONTINUF 

326 

C 

327 

C 

♦*♦•»  CALCULATION  OF  OPEN-LOOP  EIGENVALUES  ♦»»»♦ 

32  8 

C 

- 

329 

CALL  EIGENIN,A,24,22,DUN,0) 

330 

PRINT’I/T2,9X,A//14X,A,7X,A/)’,*0PEN  LOOP  EIGENVALUES*, 

331 

:  *RE AL*  ,* IMAGINARY* 

332 

00  219  I » l , N 

333 

PR  I  NT *11  OX, 2IE12. 5, 2X1)*, 24(1), 2211) 

334 

219 

CONTINUF 

335 

C 

336 

c 

♦♦♦♦♦  CALCULATION  OF  CLOSED/LOOP  SYSTEM  EIGENVALUES  •••• 

337 

c 

338 

CALL  n*uH8,L,N,*I,N,  RESULT) 

3  39 

CALL  «5URIA, RESULT, N,N, 21) 

340 

c 

341 

c 

*****  NOTE:  21  -  Al ONE )-B IONE )*U ZERO) 

342 

c 
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353  CALL  eICEN(N,71,/5,Z2,0UN,01 

355  PR1NT*(/T2«9X,A//15X, A, 7X,A/1*, 'CLOSED  LOOP  EIGENVALUES*. 

355  :  ‘REAL* ••IMAGINARY* 

356  00  220  1*1, N 

35  7  PRINT* ( 10X.2(E12.5.2XI  I *, 25  { 1 1 ,22( II 

358  220  CON  T  I  HUE 

359  C 

3 50  C  CALCULATION  OF  PHI,  DREV,  AND  GANNA 

351  C 

352  TAU-TINCRM 

353  CALL  *EXP(N,A.TINCR.?HI I 

355  CALL  NEXP1N.A .TAU.DREV) 

355  CALL  INTEATIN, A, RESULT, TINCRI 

356  CALL  nnuL(RESULT,B,N,N,.»,GANN*1 

357  C 

358  C  •«•*•  CALCULATION  OF  THE  ESTIMATOR  GAIN  NATRIX  ISI  *•••• 

359  C 

360  CALL  TPANSP(A.ATRANS.N.N) 

361  CALL  TRANSPIC.CTRANS.VN.N) 

362  221  CALL  I NP« TX I HANOLE , V , YN. T N , RE  SUL T , 3. NODE ) 

363  CALL  G*»tNV(YN.  YN, RESULT,  VINV.'R, 01 

365  IF  (5P.NE.YNI  THEN 

365  PRINTM/T2.A/)*,  *V  NOT  FULL  RANK.* 

366  IF  (NaOE.CO.il  THEN 

367  GO  TO  1 

368  ENOIF 

369  JAR  OR  T* l 

370  GO  TO  250 

371  ENOIF 

372  CALL  nkuL(CTRANS,VINV,N,YN,YN,RESULTI  i 

373  CALL  N"UL(RESULT,C,N,YN,N,RESULT2> 

375  CALL  -PICIN,ATRANS,RESULT2.M,P,21,IERI 

375  IF  (1ER.EQ.il  THEN 

376  IF  (NOOE.EQ.il  THEN 

377  GO  TO  1 

378  ENOIF 

379  JARORT*! 

380  GO  TO  250 

381  ENOIF 

382  CALL  NNULIP, RESULT, N.N.YN, SI 

383  C 

385  C  *****  CALCULATION  OF  PSI  *•••• 

385  C 

386  CALL  INTE AT(N, A, RESULT, TINCRI 

387  CALL  N"UL(RESULT,S,N,N,YN,PSII 

388  C 

389  C  ♦♦♦••  CALCULATION  OF  THE  COVARIENCE  RA  TR ICES  **••• 

390  C  OF  THE  STATES  AND  THE  NE A  SURE  IE  NT S. 

391  C 

392  CALL  Nnu(.<0OEV,P,N,N,N,RESULTI 

393  CALL  NEXP(N,ATRANS,TAU,2ll 

395  CALL  N«UL(RESULT,Zl  ,N  ,  N,  N  ,  R  E  SUL  T3  1 

395  CALL  INPNTXIHANOLE ,  N  ,  N  ,N  .  R  E  SUL  T , 3.N00E I 

396  CALL  INTECIN, A, RESULT, RESULT2.TAU1 

397  CALL  NAOO(RESULT3,RESULT2,N,N, RESULT! 

398  CALL  ■"*UL(0REV,P,N,N,N,RCSULT2I 

399  CALL  nnUL(RESULT2.CTRANS,N,N,YN,RESULT3I 
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*00  C4LL  NM<jl<RESULT3.V!NV,N,YN,VN,RESULT2  > 

*01  CALL  MmuHRESULT2,C,N,YN,N,RESULT31 

*02  CAUL  M~UL(RESULT3,P,N,N,N,RESULT2> 

*03  CALL  NMULIRESULT2.Z 1, N , N , N , RESUL T 3  I 

*0*  IOL-.001 

*05  CALL  TRANSP<2,RESULT2,N,NI 

*06  CALL  MLINE0<N, RE SULT2, RESULTS, 21, T0L1 

*07  CALL  *A00<RESULT,Z1  ,N  ,  N ,  R  ES  UL  T  3  > 

*08  CALL  m*uL(C.RESULT3,YN»N,N, RESULT! 

*0*2  CALL  MMUL(RESULT,CTRANS.YN.N,YN,RESULT2 I 

*10  C 

*11  C  «••♦*  NOTE :  2  -  A( BAR  )  -  A  1 2ER 0 > -8 1 2 E R 0 >>L I  2 E RO >  *••*♦ 

*12  C  ««**•  RESULT2  -  DISPLAY  VARIABLE  COVARIANCE  MATRIX  **«•• 

*13  C  *♦♦«*  RESULTS  -  STATE  COVARIANCE  MATRIX  ***** 

*1*  C 

*15  C 

*16  C  «*«**  COMPARING  GUESSED  NOISES  TO  VAR  I ENCE  S 

*17  C 

*18  GUESS-0 

*10  00  223  I -  1 , YN 

*20  SUMl-S0RT(R£SULT2<l.I I *2 » 

*21  SU-l-THRSHOI I 1/SUM1 

*22  SUM2-ERFCI SUMil 

*23  SUMl-SUM2**2  * 

*2*  RESULT | l,I)-.031*15O*(RESULT2I I »! 1/SUN1I 

*25  SUM1-A8S((V(1,I}-RESULT(I,I1  I  7  RE SULT (  I  ,  I  )  I 

*26  IF  ISU'«I.LT..05>  THEN 

*27  GO  TO  223 

*28  ENOIF 

*20  VU.II-RESULTt  !,!l  v 

*30  GUESS- 1 

*31  22 3  CONTINUE 

*32  00  22*  tMN-NN)~IN-U,(N-NM! 

*3  3  211  1  ,1  1  » .0  00*2*8*RESUL T 3 ( I , I  I 

*3*  VI I  ,!  >*WI  I  ,I)*(TAUNU  ,  I  l**2» 

*  35  SUMl - A8S I  I M It , 1 1 -2 1 II , I > ) 72 1 1  I , 1)1 

*36  IF  <SUMl.LT,. 051  THEN 

*37  Ml  l  ,  I  I-Wll,  Il/C  TAUNC  I  ,  I )  **2  ) 

*38  GO  TO  22* 

*30  ENOIF 

**0  Nil, tl*Zl< I, II 

**l  Mil, !)-MII,I)/ITAUNlI, 11**21 

M*2  GUESS-l 

**3  22*  CONTINUE 

***  C 

**5  C  *****  note:  RESULT  -  I. 01  P I *0 l SPLAY  CO V AR I ANCE » 7  I F 1  HA T  )) »*2  **♦•• 

* *6  C  •••••  21  -  .003  PI*STATE  COVARIANCE  •••♦. 

**  7  C 

**8  IF  I  CUE  5S.E0.il  THEN 

**0  CO  TO  221 

*50  ENOIF 

*51  C 

*52  C 

*53  PRINT*, 'STATE  COVARIANCE  MATRIX* 

*5*  DO  230  I - 1 , N 

*55  PRINT* <10I3X,E10.*I1’*1RESULT3II,J)»J-1»N1 

*56  230  CONTINUE 
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♦  5  7  250  PS INT* I /T2.A/1  • .  *ThE  P  MATRIX  IS:' 

45*  DO  241  I-l.N 

459  PRINT*! 101 3X.E10.4II *, (PI  1,41 , J-I.N) 

460  241  CONTINUE 

461  PRINT  *I/T2,A/|*,  •  THE  KALMAN  FILTER  CAIN  MATRIX  (S3  IS:  • 

462  00  242  I-l.N 

46  3  PRINT'dOl  3K.E10.4I  I*.  (S(  I.  Jl  .  J-l.TNI 

464  242  CONTINUE 

465  C 

466  C  *****  CALCULATION  OF  CLOSEO-LOOP  ESTIMATOR  EIGENVALUES  *•♦*» 

467  C 

468  CALL  MMUL(S.C.N,rN,N,2I 

469  CALL  MSU8<A,2.N.N,21| 

470  C 

471  C  *****  NOTE:  21  -  A10NEI-S*C 

4  72  C 

473  CALL  EICENIN. 21.24.22. DUM.OI 

474  PRINT*t/T2,9X.A//14X,A,7X,A7)*, 

475  :  ‘CLOSEO-LOOP  ESTIMATOR  E I GENVALUES * . *  RE AL * . * l "AGI NAR Y • 

476  00  243  I-l.N 

477  PRINT* I 10X.2IE12.5.2X) I • ,24| |) ,22( (I 

478  243  CONTINUE 

4  79  C 

480  C  REFORMING  AO  ANO  BO  FOR  N?XT  RUN  *♦*♦* 

481  C 

482  III-l 

483  00  246  1 -IN-MM)— IM— 1 1 , (N— MM) 

484  00  244  I  1-1,8 

485  BII.llI-O. 

486  244  CONTINUE  1 

487  DO  245  J.|N-MM)-(M-1 >, (N-MMI 

488  AII.JI-O. 

489  245  CONTINUE 

'490  SII.IlII-l. 

491  III-III-I 

492  246  CONTINUE 

493  C 

494  C  THEN  RETURN  TO  START  FOR  THE  NEXT  PROBLEM.  ♦**•* 

495  C  ON  INTERCOM,  R£AO  REMAINING  CAROS  FOR  THIS  PROBLEM, 

496  C  »••••  IF  ON  INTERCOM  ANO  JABORT-1,  RETURN  TO  START.  IF  NOT  *»••• 

497  C 

498  250  IF  (KTRAJ.E0.21  THEN 

499  GO  TO  260 

500  ELSE  IF  (KTRAJ.EO.OI  THEN 

501  GO  TO  290 

502  END  IF 

503  C 

504  C  PLOT/LIST  INFORMATION  »••»* 

505  C 

506  IF  (MODE. EO. II  THEN 

507  PRI NT*< /T2.A/T2, A, 12, A/ )• ,  'ENTER  TITLE  OF  PLOT,*, 

508  :  •(•,30,*  CHARACTERS  MAX  I . • 

509  ENOIF 

510  READ  (5.*<3A10|*>  TITLE 

511  260  IF  (M00E.E0.1I  THEN 

512  PRINT* (/T2 ,A/TS, A/T8.2 kf I* , 

513  :  'ENTER  THE  FOLLOWING  INFORMATION  FOR  PLOT/LIST:*, 
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SI* 

'PLOT  TI*E  INCREMENT  (DT)  AS  A  REAL  NUHBER 

5  IS 

•NU»8£R  OF  OUTPUT/STATE  T I HE  RESPONSES  DESIRED 

516 

'FOR  EACH  PLOT,  (HAX  IS  FOUR).' 

517 

END  IF 

515 

READ*,  OT.NOUT 

519 

IF  (NOUT.OT.*!  THEN 

520 

PRINT' ( /T2, A/1 WAKE  UP,  EIGHT-BALL.  I  SAID  THE  HAX  IS  FOUR. 

521 

IF  MOOE.EQ.l)  THEN 

522 

GO  TO  260 

523 

ENOIF 

52* 

GO  TO  265 

525 

ENOIF 

526 

IF  (NOOE.EO.l)  THEN 

527 

PRINTM/T2,2A,I1,A/T2,A,U,2A2T2,2A.II,A/T2,2A)', 

525 

•FOR  EACH  STATE  T I *E  RESPONSE  DESIRED,  ENTER 

529 

•A  • , l , '  X  N  HATPIX  (COUT)  OF  2ER0S  WITH', 

530 

•A  "•.l.*"  PLACED  IN  THE  POSITION  CORRESPONDING 

531 

•TO  THE  DESIRED  STATE.* , 

532 

•FOR  EACH  OUTPUT  TINE  RESPONSE  OESIREO  ENTER 

533 

•A  ',1,*  X  N  HATRIX  (COUT)  CONSISTING  OF  THE 

53* 

•ROW  OF  THE  C  HAT R I X  CORRESPONDING  TO  THE  *. 

535 

•OESIREO  OUTPUT  1 Y  > . • 

536 

ENOIF 

537 

2  65 

00  266  1*1, NOUT  * 

538 

REAO  15, ♦!  ICOUT I I , J ) , J*1 , N  ) 

539 

2  66 

CONTINUE 

5*0 

PRINTM/T2.A/)',  'THE  PLOT/LIST  OUTPUT  HATRIX  ICOUT! 

IS:  * 

5*1 

00  267  I » l , NOUT 

5*2 

PRINT* 1 101 3X,610.*n* , ICOUT! I,J!,J*1,N) 

5*3 

267 

CONTINUE  ‘ 

5** 

IF  IKTRAJ.E0.2I  THEN 

5*5 

GO  TO  280 

5*6 

ENOIF 

5*7 

IF  IHOOE.EO.l)  THEN 

5*8 

PRINTM/T2.A/1*,  'ENTER  LABELS  FOR  EACH  OUTPUT.' 

5*9 

ENOIF 

550 

268 

00  275  1*1, NOUT 

551 

GO  TO  (270,271,272,273), I 

552 

270 

REAO  1 5, • I  2 AiO) • 1  Y IT  I TLE 1 1 > , Y1 T I TL E I  2 ) 

553 

GO  TO  275 

55* 

271 

REAO  (5,*I2A101'I  Y2TITLEIH,  Y2TITLEI2) 

555 

CO  TO  275 

556 

272 

REAO  I5.,(2A10!‘I  Y3TITLEI1)  •  Y3TITLEI2) 

- 

557 

GO  TO  275 

558 

273 

REAO  (5,*(2A10I*I  Y*TITLE(1).  Y*TITLE(2I 

559 

275 

CONTINUE 

560 

C 

561 

C 

562 

280 

NPOINTS-ITF— TOI/OT 

563 

NLEN-300 

56* 

IF  (NP0INTS.GT.298. )  THEN 

565 

NPO I  NTS*  29  8 , 

566 

ENOIF 

567 

I  I -0 

568 

COUNT-O 

569 

L  IHIT-OT/TINCR 

* 

570 

290 

IF  (JABORT.EO.il  THEN 
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571 

CO  TO  802 

572 

ENOIF 

573 

C 

5  75 

C 

*  ♦ 

575 

c 

♦ 

DEFINING  OPTI-al  PILOT  V  AP  I  ABL  E  S 

576 

c 

* 

UOPT: 

PRESENT  0PTI8AL  INPUT 

577 

c 

# 

UOLD: 

PAST  OPT!«AU  IMPuT 

578 

c 

• 

UOELAY: 

INPUT  DELAYING  MATRIX 

579 

c 

• 

OUOPT: 

OELAYEO  OPTI-AL  input 

580 

c 

♦ 

XNf  w: 

PRESENT  STATES 

581 

c 

# 

XCLO: 

PAST  STATES 

582 

c 

• 

rwew: 

PRESENT  DISPLAY  VARIABLES 

583 

c 

♦ 

YOELAV: 

0 I S  PL  AY  VARIABLES  DELAYING  MATRIX 

588 

c 

♦ 

Y  : 

delayed  display  variables 

585 

c 

* 

TEST: 

ESTIMATOR  DISPLAY  variable  error 

586 

c 

* 

DxhaTn: 

PRESENT  ESTIMATED  OELAYEO  STATES 

587 

c 

• 

DXHATO : 

PAST  ESTIMATED  delayed  STATES 

588 

c 

* 

0 XH  A  T  N ; 

PRESENT  PREDICTED  STATES  DUE  TO  OPTIMAL 

589 

c 

• 

INPUTS  (PREDICTED  OPTIMAL  STATES) 

590 

c 

♦ 

OXHATO: 

past  predicted  OPTIMAL  STATES 

591 

c 

♦ 

XDELAY : 

PREDICTED  OPTIMAL  STATES  DELAYINC  MATRIX 

592 

c 

• 

OOXHAT : 

OELAYEO  PREDICTED  OPTIMAL  STATES 

5<>3 

c 

• 

OEE  : 

DELAYED  STATE-  Rf  CONSTRUCT  I  ON  ERROR 

598 

c 

• 

EE: 

STATE  RECONSTRUCTION  ERROR 

595 

c 

• 

xhat  : 

RECONSTRUCTED  STATES 

596 

c 

♦ 

P  F  * : 

PERFORMANCE  INDEX  (ONE  ITERATION) 

597 

c 

♦ 

PF  NT: 

TOTAL  PERFORMANCE  INDEX 

598 

c 

♦  • 

*99 

c 

600 

c 

V 

601 

c 

♦  ♦ 

INITIALING  OPTICAL  PILOT  VARIABLES 

602 

c 

60  3 

PF--0.0 

608 

pfmt-o.o 

605 

*3 

-0. 

606 

U3 

*0. 

607 

00 

292  I • l • NO  I 8 

60  8 

00  291  J-l.NOIM 

609 

X2I I , J) *0 

• 

610 

U2  1  I . J) -0 

• 

611 

YOELAYC  I  , 

Jl-O. 

612 

XOELAYI  I  , 

Jl-O. 

613 

UOELAYI  I  , 

Jl-O. 

618 

291 

CONT INUE 

615 

292 

CONTINUE 

616 

00 

297  I  ■  1 .  N 

617 

YESTIM-O. 

618 

UOLD I I »-0. 

619 

UOPTI I 1*0. 

620 

XNFVM  ll-O. 

621 

OXHATNI I ) *0. 

622 

OXMATOI I 1*0. 

623 

OFE  1  I  1  *0. 

628 

EE  I l 1*0. 

625 

XHAT ( I 1*0. 

626 

OXHATNI I|*0. 

627 

XOLOU  1-0. 
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628 

oxha  ro( 1 1 -o. 

629 

OUOP  T ( ll-o. 

630 

OOXHAT ( I ) -0. 

631 

YNEw ( 1) -0. 

632 

Y( I 1 -0. 

633 

297 

CONTINUE 

63* 

T0L0-T1 

635 

CALL  INP«TXIHAN0L£.XI.N,1,X0L0.3.,,00E> 

636 

CALL  INP"TXIHANOLE,XI,N,1,DXHATO,3,NODE) 

637 

CALL  INP8TX I HANOLE  ,  Ut  ,8,1 ,UOLO, 3,-ODE ) 

638 

C 

639 

C 

*****  THE  OPTIMAL  pilot  loop  ***** 

6*0 

C 

6*1 

300 

TNE'A»TQLO*T  INCR 

6*2 

IF  ITNEW.GE.TF)  THEN 

6*3 

CO  TO  *00 

6** 

ENOIF 

6*5 

C 

6*6 

c 

««♦*♦  CALCULATION  OF  THE  OPTIMAL  INPUT 

6*7 

c 

6*8 

CALL  SmmulTIL.-I. , " , N , RE  SUL T ) 

6*9 

CALL  88ULIRESULT.XHAT.8,N,1 .UOPT) 

650 

c 

_ 

651 

c 

*****  DELAYING  THE  OPTICAL  INPUT  BY  TAU 

SECONOS  ••• 

652 

c 

653 

00  311  J-1,8 

65* 

00  310  1-1,8 

655 

UOELAYI  I,J)-U0ELAY(I,J-1> 

656 

310 

continue 

657 

311 

CONTINUE 

V 

658 

00  312  1-1.8 

659 

UOELAYI I ,91 -UOPT I I 1 

660 

312 

CONTINUE 

661 

00  313  1*1.8 

662 

OUOPTI I ) -UOELAYI 1,1) 

663 

313 

CONTINUE 

66* 

C 

665 

C 

CALCULATION  OF  SYSTEM  STATES  ♦♦♦•* 

666 

C 

667 

CALL  "MULIPHI  .XOLO.N.N, 1, RESULT) 

668 

CALL  8  8UHCA8HA,U0PT,N,8,l,RESULT2> 

669 

CALL  MAOOIRESULT.RE  SULT2 , N , 1 , XNE8 > 

670 

c 

- 

671 

c 

♦♦«««  CALCULATION  ANO  DELAY  OF  DISPLAYED 

VARIABLES  1 

672 

c 

673 

00  121  J-1,8 

67* 

00  320  I-l.YN 

675 

YOELAYI I.JI-YOELAYI I.J-l) 

676 

320 

CONTINUE 

677 

321 

CONTINUE 

678 

CALL  5«UL(C,XNE8,YN,N,l,YNEM) 

679 

00  322  I-l.YN 

680 

YOELAYI I ,91 -YNEWI l ! 

681 

322 

CONTINUE 

682 

DO  323  I-l.YN 

683 

YIII-YOELAYIl.l) 

68* 

32  3 

CONTINUE 
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6*S  C 

6*6  C  ♦♦♦•♦  CALCULATION  Of  ESTIMATED  STATES 

6*7  C 

688  CALL  -*UL  (C  .CXHATN,  YN  ,N,  1  ,  RESUL  T  I 

68  9  DO  110  I-l.VN 

600  VEST  I  I )-Y( 11-RESuLTI l , l> 

691  310  CONTINUE 

692  C 

691  call  ■*UL(PHI .OXHATO.N.N, 1, RESULT) 

69*  CALL  -*UHGA-><a,DU0PT,n,m,1,RESULT2> 

695  CALL  "AOOIRfSULT.RESULT^.N, 1, RESULT) 

696  CALL  ■"‘•ULCPSI, TEST, N.YN.l, RESULT?) 

697  CALL  *AOOI RESULT .RESULT2 ,N, l .OXHaTN) 

69*  C 

699  C  ••*«*  CALCULATION  OP  PREDICTED  OPTICAL  STATES  *•••♦ 

TOO  C 

701  CALL  N*UL  (PH!  .OXHATO  .N  ,N,  1  .RESUL  Tl 

702  CALL  m*UL(GA*1*A, uOP  T ,  N  ,  N , 1 .RESULT2 ) 

701  CALL  ■'A00(RESULT,RESULT2,N,1,0XHATN| 

70*  C 

705  C  ♦  ♦♦•«  OfLAT!NG  PREDICTED  optimal  states  ***** 

706  C 

707  DO  3*1  J  "  1 , 8 

708  00  1*0  I-l.N  • 

709  XOfLAYl I ,J)-XDELATII , J*l) 

710  3*0  CONTINUE 

711  3*1  CONTINUE 

71?  00  3*2  1*1. N 

711  XOELAYI l  ,9)-0xHATN( I) 

71*  3*2  CONTINUE  «. 

71 5  00  1*1  I-l.N 

716  OOXHAT ( I ) -XOELAYI I , l ) 

717  3*3  CONTINUE 

71*  C 

719  C  ••♦*•  CALCULATION  Of  RECONSTRUCTED  SYSTEM  STATES  (XHATI 

720  C 

721  CALL  "SUPIOXhaTN, OOXHAT, N.l, DEE) 

722  CALL  --ULIOREV.CEE.N.N.l.EE) 

721  CALL  MAOOIEE.OXHATN.N.l.XHAT) 

72*  C 

72 5  C  ♦«•♦*  UPOATING  VARIABLES  FOR  NEXT  ITERATION  *••♦• 

726  C 

727  TOL 0 • TNC  W 

728  00  3  SO  I-l.N 

729  XOLOf  I  l-XNEWI  I) 

730  OXHATOI I  1 -OXHATNI I ) 

731  OXHATOI I ) -OXHATnI I ) 

732  3  SO  CONTINUE 

713  00  1SI  I-l.N 

71*  UOLO (  I  1 -UOP  T ( I  1 

73S  3SI  CONTINUE 

716  C 

717  C 

71*  C  ••••♦  CALCULATION  of  THE  PERFORMANCE  MEASURE  ***** 

719  C 

7*0  CALL  TRAnSPI XCLD.RE SULT.N.l 1 

7*1  CALL  NMULIXOLO, RESULT, N.1.N.X2) 


GAE/AA/SOD-3 


742 

CALL  08«UL( X2,0,N,X3) 

743 

CALL  TRAN'S?  ( UOLO.RE  5ULT2. 4,  1 1 

744 

CALL  M MUL (UOLO.RE SUL T 2, M, l ,8,u2 ) 

744 

CALL  0«MUL(U2  ,R,8,U3I 

7*6 

PF  9 *X  3  *U  3 

747 

P  p  «  T  .  p  f  i  T  ♦  P  F  N 

748 

C 

749 

IF  (KTRAJ.EO.O)  THEN 

740 

GO  TO  300 

741 

ENOIF 

742 

c 

743 

c 

. 

74  4 

c 

STORAGE  OF  1 NF OR 8A T I  ON  F OR  PLOT/LISTING  •♦••• 

744 

c 

746 

COUNT-COUNT-l 

747 

If  (COUNT. EG. 1)  THEN 

748 

I  I  -  I  1  ♦  1 

74** 

TP ( I I 1 -TOLD 

760 

CALL  8*UL 1 COUT.X OLD, NOUT ,N, 1 , OUTPUT) 

761 

i noe  x *  1 1 

762 

DO  360  1*1, NOU T 

76  3 

YA( INOE  X l-OUTPUT ( I 1 

76* 

INDEX-INO£X«NLEN 

764 

360 

CONTINUE  * 

766 

ENOIF 

767 

IF  (COUNT. CE.LI8ITI  THEN 

768 

COUNT -0 

769 

ENOIF 

770 

c 

771 

c 

*«♦♦♦  end  OF  OPTIMAL  PILOT  LOOP  •♦•••  * 

772 

c 

773 

GO  TO  300 

774 

c 

774 

c 

♦  ♦♦••  PLOT/LISTING 

776 

c 

777 

4  00 

PRINT*  ( /T2.A.E20.6/ 1  * .'THE  VALUE  OF  THE  PERFORMANCE  INDEX  IS:  *, 

778 

:  PF  *  T 

779 

402 

IF  (KTPAJ*FO.O)  THEN 

780 

GO  TO  420 

781 

ENOIF 

782 

IF  (JA80RT.E0.il  THEN 

783 

GO  TO  410 

784 

ENOIF 

784 

IF  (KTRAJ.EO. l.OR.KTRAJ.EO.  3)  THEN 

786 

CALL  DR  AW  2  (  TP  ,YA,YA(  301)  ,YA(  6011  ,YA  (901  1  ,  NOUT  .  NOUT  .  NP  0 1  NT  S  ,  0, 1) 

787 

ENOIF 

788 

IF  (KTRAJ.EO. 2. OR. KTRAJ.EO. 3)  THEN 

789 

IF  ("OOE.EO.ll  THEN 

790 

IG*(NP0InT5*2)/3 

791 

ELSE 

792 

IG-(NP0INTS*41/5 

793 

ENOIF 

79* 

00  404  LA-1, NOUT 

794 

IF  (M00E.E0.1I  THEN 

796 

PR INT*( /T4,3( A, 10X.A, 11 ,9X1/1 • ,*T* ,* Y*,LA 

797 

ELSE 

798 

PRINT*(/T4,4IA,10X,A,Il,9X)/)*,*T,,*Y*,LA 
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799 

enoip 

800 

[NC-(L4-1 1  *NLEN 

80 1 

oo  403  i e- 1 ,  i c 

802 

PB|NTM12,5(P8.3.F12.6.6Xll,t 

803 

;  (  TP{  10  1  .TA(  IOMNCI  .  10  -  IE.NPOINTS  ,  Id 

SO* 

603 

CONTI  HUE 

805 

<»0<. 

CONTINUE  k 

806 

E  NO  I E  - 

807 

C 

808 

c 

*•♦•*  IP  NOPE  PLOTS/L  I S  T I NC  S  POP  THIS  PR08LE8  •  *♦•♦ 

809 

c 

•  •••♦  APE  0E5IPE0.  3  E  T  ""0PE"-2.  •  •»*• 

810 

c 

811 

610 

IF  (800E.E0.1I  THEN 

812 

PP  INTM/T2  ,A/T2,  A,  I  1/1  •  .  MF  POPE  PL OT S/ L I S T I NCS  FOP*. 

813 

:  'THIS  P  P  08  L  £  6  APE  OESIBEO,  ENTEP  A  ',2 

81*. 

ENO  IF 

815 

PE  AO • .  "OPE 

816 

IF  ( JA80PT. EO. 1 )  THEN 

817 

CO  TO  620 

818 

END  I  F 

816 

IF  ("0PE.E0.2I  THEN 

820 

CO  TO  250 

821 

ENO  I  F 

822 

C 

823 

C 

♦••••  IF  CHANCES  TO  THE  STS TE"/PP08LE 8  APE  TO  BE  *AOE  AND 

826 

c 

*****  a  NEW  PP08LE*  CPEATEO.  SET  "CHANCE"-!.  ♦•••' 

825 

c 

826 

6  20 

IF  ("OOE.EO.ll  THEN 

82  7 

PPINTM/T2,A/T2.A,!1/I  *  ,  MF  CHANCES  TO  THE  SrS  •  "/ PP  OBlE  * 

828 

t  *APE  TO  BE  8A0E  AND  A  NEW  PR08LE "  CREATED.  ENTER  A  •, 

826 

ENO  1  F 

830 

READ*. CHANCE 

831 

IF  (CHANCE. €0. 11  Then 

832 

CO  TO  l 

833 

ENO  I  F 

83* 

C 

835 

500 
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P=00T31H, A(I> ,S) 

22 

L  =  I 

23 

00  40  J=1,M 

24 

A(L)=A(L)-0(J)»P 

25 

40 

L*L*N  CIM 

26 

f* 

«# 

COLUMN  MODIFICATION  OF  ' 

27 

IF(NR  .EQ. u)  GO  TO  63 

28 

KJ=RI *N-i 

29 

00  55  I=KI,KJ 

30 

P  =  OOT  3 (M, 0 ( I ) , S) 

31 

L  =  I 

32 

00  50  J=i,H 

33 

0(D=  r.(L)  -0(J)  »F 

34 

5  C 

L=L*NOIM 

35 

60 

IF(<a  .LE.KK)  RETURN 

36 

0(1)  =  -T 

37 

00  73  I*2,M 

38 

7  C 

0 ( I)  =  9, 

39 

RET'JR  W 

43 

END 
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SUBROUTINE  EIGEN 


1 

SUBROUTINE  EIGENfN.AA.RR.RI.O.NR) 

z 

01  PENSION  AA(1),RR(1>  ,RI  (1),  0(1)  ,  0R(3> 

3 

CONNON/HAINI/NPIM.HOIM,  A  <  1)  'INO-l'KIN,  KCUl 

4 

E1UI VALENCE  (P,QR<1)  )  ,<B,OR(2)>,(i,OR<3>> 

5 

GAMA} al.E-7 

6 

NN  =  N*  <N-1)*NDIM 

7 

DO  2  1=1, N 

8 

00  1  0=1,  NN,  NOIN 

9 

1 

A  <  J)  =  AA ( J) 

10 

2 

rich  =o.r 

11 

IF(N. LE.2)  GO  TO  30 

12 

NE  =  ( N  -2)  »NDIN 

13 

H=N-1 

14 

KI=MOIN 

15 

KN  =  M 

16 

00  20  KK=  2  ,  NE,  HOI  N 

17 

KL  =  1*  KK 

18 

KM=<N*NOIM 

19 

X=ASS  (A  OOO  ) 

20 

Y*X»G  A«AO 

21 

00  5  I=KL,KN 

22 

5 

X  =  X*ABS(A(I)J 

23 

IF(X. GE.T)  GO  TO  15 

24 

DO  5  1=KL,KN 

25 

6 

A(I»  =  ".0 

26 

GO  TO  16 

27 

15 

CALL  KEOCKN.A.KK.KI.KN,  N.AOOO  ,NR,  01 

28 

16 

KN  =  KM 

29 

KI  =  KI  4N0IH 

30 

2  u 

M=i-i 

31 

C - 

A  IS  IN  UPPER  HESSEN 9ERS 

32 

3  C 

MN  =  NN 

33 

M=N 

34 

TT  =  0. 

35 

NE  =  NN  -NO I M 

36 

C 

NEXT  ITERATION 

37 

4C 

NA=NE-HOIM 

38 

ns=nm-moim 

39 

gana=  c-amao 

40 

45 

ITS  =  fl 

41 

X  =  MM  ►> 

42 

IF(N. EC . 1 )  GO  TO  210 

43 

C 

LOCK  TOR  ZERO  ON  SUB  -  0 IAGONA 

44 

50 

J-M1 

45 

SO 

TOL  =  GAMA»  (  ABS  (  A  (  J  )  )♦  ABS  (  A  (J-NOID  ( 

46 

IP(A3S(AtJ-N0IH) ) .LE.TOL)  GO  TO  7! 

47 

J= J-HOIH 

48 

Ie(J.GT.l>  GO  TO  60 

49 

7  C 

X  =  A  ( H  V ) 

50 

Ic  ( J. EO.NM)  GO  TO  210 

51 

Y  =  A  ( N  S) 

52 

W*  A  (  N  E)»A  (HN-1) 

53 

IR(J. EO.NS)  GO  TO  220 

54 

I'-(ITS.EO.SO)  GO  TO  300 

55 

IFUTS.NE.20.ANO.  ITS.NE.1V>  GO  TO  )) 

56 

TTUT  »X 

57 

00  81  1=1 , HH, HOI M 
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SUBROUTINE  EIGEN 


58 

59 
6C 
61 
62 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 
89 
99 

91 

92 

93 

94 

95 

96 

97 

98 

99 
130 
101 
102 

103 

104 
1  05 
1C6 
107 
1  08 
1  09 
110 
111 
112 
113 
1  14 


80  »(n=MI)-X 

Z  =  A3S  ( A  t  N£ )  )  *A 9 S  (  A  (NA> )  , 

x=.z ;•? 

v=x 

Wr-,4  375*  7*  2 
9C  I  TS=  1 *ITS 

C  10  <K  FOR  TWO  CONSECUTIVE  SMALL  5U9-OIAGONAL  ELEMENTS 

KiNS-  ►DIM 
IS C  7  =  A(K  ) 

RiX-’ 

S  =  r-’ 

P=<R‘S-W)/AT1*K>*ATK*NDIH> 

8=  A  ( *  ♦MOIH) -Z-4-S 
RiA<<  *NOIM*2) 

IF (<. EO. J!  GO  TO  115 

TOL  =  3  AM&»A8S(F>»  TABS ( A< K-M3I M) ) » A e> < Z ) ♦ A  9 S ( A  « *M 01  M) )) 
IF  (A9S  (A  (K-NOIN)  )  «  <A9S(  B)  *ABS  (R)  )  .  .E.  TOL>  GO  TC  113 
K=K-MOIN 
GO  TO  ICC 

C  9EGIN  ITERATION  WITH  DIAGONAL  ELEMENT  K 

ltC  KI=<<»NOIM*1>/MOIM 
L=3 

KL=M-1*KI 

DO  230  KK=K,NS,MOIM 
IF(«.EQ.NS>  L  =  2 
KC  =  «  *L 

IF  I<C .G  T . KL )  KC- KL 
ICT<K  ,N£.  K)  GO  TO  150 
IF(J.KE.K)  ATK-NOIM)  i-A  (K-NOIM) 
call  5E0CT<N,Af<K,KI,KC»Lf9R,  NR,  31 
GO  TO  175 
15  0  KN  =  << -NOI M 

cal.  feo;ttn,a,kn,ki  ,k:,l,at<n>  ,nr,3> 

175  KL=KL*N0IM 
2CS  KI  =  KI*NDIM 
GO  TO  50 

3  c 0 UNO  ONE  ROOT 

210  A  (  MM)  iX»TT 
RR  <  Ml  iX»TT 

IF(1.  E0.1)  GO  TO  323 

HiM-1 

A(NE)  =9, 

NE  =  NA 
MM  sNS 
GO  TO  40 

0  FOUND  TWO  ROOTS 

2  23  P-.  5*  <T-X) 

RiO»B  <W 

7  =  S0R  T  TABS (9)1 
XiX»T  T 
A  <MM)  iX 
A  (  Nil  =Y  ♦  T  T 

IF  T9.  LT.f.)  GO  TO  29  0 
:  RE  «L  ROOTS 

TiR*SIGn<7,P> 

241  RR(M-l)iX«x 
RR ( M I  =X - w / T 
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S URROUT [ME  EIGEN 


115 

X=4(N  E> 

116 

* 

RECUCE  TC  DIAGONAL 

117 

Y=S0RT(x*X*7»Z> 

118 

D 

RE  CUCE  TO  DIAGONAL 

119 

Y  =  SDRT  <X*X»Z*Z) 

1  ZD 

PsX/Y 

121 

9-7/Y 

122 

2 

RD  L  AND  COL  JMN  HOD 

123 

DO  250  I=NSiNNiNOIM 

12*. 

7=4(1  ) 

1  25 

Y  =  A  (I  .1) 

126 

4(I)*9*Z»P»Y 

127 

25  3 

A  ( I  ♦  1  )*8*Y-P*Z 

128 

<C  =  NE -H*l 

129 

DO  260  I=KC,NE 

130 

7=4(1) 

131 

Y  =  A(I  ♦NOIM) 

132 

All)  =  E»Z*P»Y 

133 

263 

A(I»NCIM) =9»Y-P»2 

134 

* 

DOLUMN  4  OOI PIC  ATI  ON 

135 

IMNR  .EQ.O)  GO  TO  275 

136 

«  =  <D  *N-1 

1  J7 

00  27  t  I=KC,KK 

1  38 

7*0(1) 

139 

Y=0(I  *NOIM) 

140 

0(1)  =  9*7*P»Y 

141 

275 

0(I*NDIN)=S»Y-P*7 

142 

275 

M=4-l 

143 

A(  44)  =PR(M»1) 

144 

4  <NE)  =0. 

145 

280 

:p<4.  EO. 1 )  GO  TO  320 

146 

4=4-1 

147 

A  (  N4  )  =0. 

148 

NE  =  N4  -4014 

149 

M4  =  NS  -401  4 

150 

GO  TO  40 

151 

C 

COMPLEX  ROOTS 

152 

231 

RR (4) =X»P 

153 

RI  (4)  =-Z 

154 

4=4-1 

155 

RR(4) =  X*P 

156 

RI (4)  =’ 

157 

GO  TO  280 

158 

159 

30  C 

WRITE  (XOJT.31C) 

160 

31  C 

P0R44K41H;  ONF  digit 

161 

G444=  G4M4»lt . 

162 

GO  TO  45 

163 

320 

CONTI NUE 

164 

RE  T  JR  N 

165 

END 
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SUBROUTINE  MRIC 


1 

SUBROUTINE  PRI C  (  N  ,  A,  S  ,0  ,  X  ,  7,  IER) 

2 

DINENSION  A(l) ,S(1),0(1) ,X(1) ,7(1),  T?(80).,TIPES(3>  ' 

3 

REAL  76(12)  ,  Z7  (12)  ,78  (144)  ,7  9  (12) 

4 

INTEGER  EIG 

5 

CONNO  MPA  INI/  NO  !•-  ,N0IP1  ,F  (1 )  /MAIN’/G  (1)  Y  I  NO  U/ KIN,  KOOT 

6 

CON1ON/BOH0/EIG 

7 

OATA  11H£S/2.C,a .5 ,4.0/ 

S 

NN  =  N*  KOI “ 

9 

TOL  =  l  .E-03 

10 

1  Tl  =  -.  5*ALOG(XN3RH(h,Q),.aOCl) 

11 

XRITl.LT. 1.)  T 1=  1  . 

12 

T1  =  T1  «4.»FlOAT(N)  /(l.*XNORH(N,A)  I 

13 

KEY=9 

14 

5  K£Y  =  K  EY.l 

15 

1C  CALL  EOUATE(X,S,N,N> 

16 

CALL  1NTEG(N,A,X,7,TI> 

17 

IF  (S1G.NE.1I  THEN 

18 

GO  TO  6 

19 

ENOIF 

20 

DO  11  1=1, NN 

21 

73 (I) =0. 

22 

11 

CONTI NUE 

23 

03  22  I=1,NN,N0IM1 

24 

7 9  (I)  *1. 

25 

22 

CONTI NUE 

26 

CALL  EIGEN(N,7,Z7 ,76,76, 1) 

27 

00  23  1=1, N 

28 

79(1)  =1/77(1) 

29 

23 

CONTI AUE 

30 

PRINT  •(/T2,9X,2A//14X,A,7X,A,14X,1/) •, 

31 

1  ‘EIGENVALUES  OF  CO  NT  ROLL  ABI  L  I  T  Y  /  09SERVA  BI  LI  TY  •  , 

32 

«  'GRA  MI  AN', 'REAL','  I  MAGI  NARY', 'RECIPROCAL' 

33 

00  33  1=1, N 

34 

PRINT  ,(1GX,2(E12.5,2X),10X,E12.5)',77(I),76(1),79(I) 

35 

33 

CONTI AUE 

36 

CALL  EGVEC(N,F,Z7 ,76,78  ) 

37 

PRINT  '(/T2.2A/)  'EIGENVECTORS  0-  THE  CON  T  ROLL  AB IL  ITT  /  '  , 

38 

1  'OBSERVABILITY  JRANNIAN' 

39 

00  44  1=1, N 

49 

PRINT'  (If  (3X,E  19.4))  '  ,  <7S(I»J!  ,  J  =  0  ,  NN- NOI M,  NOI P) 

41 

,  44 

CONTINUE 

42 

6 

CALL  c  ACT  OR (N, 7 , F , PR) 

43 

IF(NR.NE.N)  HRITE(KOUT,  15)  n,n,nr 

44 

15  F0RNAT(12H  GRAPNIAN  I S, I  2 , 1  NX , 1 2 , 9 N  OF  RINK  ,12) 

45 

IMID.LT.C)  GO  TO  60 

46 

CALL  CMINV(N,N,F, 7,rtR,0) 

47 

call  PAT5A(X,7,N,N,N,F) 

48 

cal.  HAT2(N,N,F,F,X) 

49 

GO  TO  18 

50 

r 

A*SX  IS  STABLE 

51 

ENTRY  PRIC1 

52 

T l =-9  9999 • 99 

53 

TOL  =  1  .E-C3 

54 

KFY=5 

55 

18  NN  =  N* ADIP 

56 

NN1=N -1 

57 

00  19  1=1, N 
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66 

67 

68 

69 

70 

71 

72 
7  3 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 
103 


00  20  J=  1 ,  1 1 
G<  J)  =  A  (J)  -f  <  J) 

2  C  J»  *  ?( Jt  *0(3) 

CALL  K.INE0<N,G,7,X,T0L1> 

cow*  r.o 

Cl -3 .  fl 

ii=i 

00  ?5  1  =  1, N 

C0N'7=  C0NV«A8S(X  (IX)-TRt  I>  > 

T  9  (I )  =X(in 

n=xi  «koimi 

25  Cl  =  :i  ♦ABS  ITRU)) 

IFCl  .GT.1.E2C)  GO  TO  53 
IF(CO 8V.GT.Cl»T0L)  GO  TO  40 
CALL  CMINV<N,N,Z,F,NR,0) 

CALL  NHUL<S,X,N,N,N,Z) 

WRITE  (KOUT , 27 ) IT 

27  F0R8A  7  ( 17  HJRICCAT  I  SOLN  IN  ,I2,1M  ITERATIONS) 

CALL  WADOKN.N.A  ,7,7,-l.a) 

IF(W»  ,  EG  •  N)RETURN 
WRITE (KOUT, 35) HR 

35  FORW&T < 26MGRICLAT  I  SOLN  IS  PS  0- - R  A  R< ,  X  3) 

RET jr N 
40  CONTINUE 
TOO  =  1  O.TOL 

WRITE  (KOUT,<-5)NAXIT,TOL 

45  F0R8a  T ( 2QM0RICCAT I  NON-CONV  IN, 12, 159  ITERS.  UP  TOL  TO  ,1°E1C.4> 
GO  TO  18 

5  C  WRITE  (K0UT,55)IT,T1 

55  F0R9AT  (29H0F  ICCAT  I  9L0W  UP  AT  I T  E  R1 TI  ON ,  I  2 ,1  2H  IINITtftL  T=,clC.5) 
IF(<ET.E0.5)  GO  TO  1 
5C  15  (<px.E0.4)  I ER= 1 
IF  KET.E  Q.4)  RETURN 
T 1  =  T i  «TIMES<KEY1 
WRITE  (TOUT  ,  6  5  )  T  1 

65  F0RiAT(14HCRESET  WITH  T=,F10.5) 

GO  TO  5 
ENO 
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THE  VALUE  OE  THE  PERFORMANCE  INOFX  t'S:  .Z9A013E*OA 
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Appendix  D 


CAL SPAN  Data 


GAE/AA/80D-3 


CHARACTERISTICS.  CONFIGURATION  2- 


GAE/AA/80D-3 


GAE/AA/80D-3 


CONFIGURATION 

. 

2-1 

WSP  ’  2-3  ?Sp  • 

FLIGHT/P1LUT 

--/  -  /" 

1883/A 

PILOT  RATING: 

OVERALL 

APPROACH  3  PIO  1 

SP  1 

W1ND/X-W1ND: 

0/04 

TURBULENCE:  None  =0.26 

FEEL : 

«  Forces;  Little  heavy,  not  uncomfortable,  okay. 


•  Displacement;  No  comments. 

•  Sensitivity;  No  comments. 


PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Initial  response  okay,  some  tendency  to  overshoot,  not  a 

problem. 

•  Predictability:  Tendency  to  overshoot  one  time  when  changing  pitch  atti¬ 

tude  aggressively.  For  most  of  the  task  it  was  not  a 
problem. 

a  Special  Inputs:  No  comments. 

•  PIO  Tendency:  No  comments. 

AIRSPEED  CONTROL:  Adequate. 


PERFORMANCE: 

•  Approach  Tasks : 

ILS:  Satisfactory. 


Visual  (Sidestep):  So  coments. 
•  Landing  tasks:  Not  done, 

a  Differences:  No  comments. 


WIND  AND  TURBULENCE:  Not  a  factor. 

SUMMARY  COMMENTS:  Annoying  overshoot  tendency  in  pitch  but  this  problem 

did  not  affect  task  performance. 


GAE/AA/Sor-3 


FEEL: 

•  Forces:  Comfortable  forces. 

•  Displacement:  No  comments. 

•  Sensitivity:  Sensitivity  okay. 

PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Responded  in  desired  fashion. 

•  Predictability:  Predictable  enough. 

•  Special  Inputs:  Very  responsive  to  turbulence,  annoying,  increased 

workload. 

t  PIO  Tendency:  No  special  inputs  or  PIO  tendency. 

AIRSPEED  CONTROL:  Pretty  good. 

PERFORMANCE: 

•  Approach  Tasks: 

ILS :  Went  well . 

Visual  (Sidestep).  Appears  okay. 

•  Landing  tasks:  Not  done. 

•  Differences:  Not  applicable. 

WINO  AND  TURBULENCE:  Turbulence  noticeable  and  is  a  factor. 

SUMMARY  COMMENTS:  Annoying  aircraft  to  fly,  hard  to  evaluate.  Tendency 

in  pitch  to  wander  off  a  problem,  tended  to  come  back. 


ll*7 


GAE/AA/SOD-3 


FEEL: 

•  Forces;  Good. 


•  Displacement:  About  right. 


•  Sensitivity;  Appropriate. 


PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Quite  prompt. 


•  Predictability:  Good. 


•  Special  Inputs:  None. 


a  PIO  Tendency:  None. 


AIRSPEED  CONTROL:  Easy,  tended  to  fly  slow. 


PERFORMANCE: 

«  Approach  Tasks: 

ILS:  Quite  good  in  smooth  air. 


Visual  (Sidestep):  Equally  good,  easy;  could  fly  more  aggressively 
than  normal. 

#  Landing  tasks:  Not  done. 


•  Differences:  Not  applicable. 


WIND  AND  TUR8ULENCE:  Smooth  generally  but  did  see  a  pitch  bucking  which  was 

undesirable  on  occasion. 


SUMMARY  COMMENTS: 


Suspect  that  it  would  be  unsatisfactory  in  moderate 
turbulence;  only  problem  was  turbulence  response. 


GAE/AA/80D-3 


FEEL: 

•  Forces:  Moderate 


a  Displacement:  No  comments. 

•  Sensitivity;  Best  compromise. 

PITCH  ATTITUDE  RESPONSE: 

•  Initial:  A  slight  delay. 

•  Predictability:  Impaired  by  mismatch  between  initial  and  subsequent 

response. 

a  Special  Inputs:  Used  ramp  type  inputs. 

a  PIO  Tendency:  None  but  a  tendency  to  overcontrol. 

AIRSPEEO  CONTROL:  Not  a  great  deal  of  difficulty;  attention  diverted 

somewhat  by  pitch  problems. 

PERFORMANCE: 

a  Approach  Tasks: 

ILS:  Worked  hardest,  reasonable. 

Visual  (Sidestep):  Reasonable, 
a  Landing  tasks:  Not  dona, 

a  Differences:  No  comments. 

VINO  AND  TURBULENCE:  Some  response  to  turbulence  but  not  a  problem;  crosswind 

—  from  right. 

SUMMARY  COMMENTS:  Major  problem  is  hesitancy  to  use  aggressive  inputs  due 

to  poor  predictability. 


GAE/AA/80D-3 


n 

CONFIGURATION 

Usp  *  2.3  fsp  -  0.57 

--/  0.5  /  - 

flight/pilot 

1889/A 

PILOT  RATING:  OVERALL  --  APPROACH  S  PIO  2  SP  J 

WIND/X-WIND:  OS/OJ  TURBULENCE:  None  M.  *  0.6a 

FEEL: 

•  Forces;  Somewhat  heavy,  but  didn't  Bind  them. 


•  Displacement;  No  comments. 

•  Sensitivity;  No  comments. 


PITCH  ATTITUDE  RESPONSE; 

•  Initial:  Okay,  no  noticeable  lag. 

•  Predictability:  Not  very  good.  Very  difficult  to  figure  aircraft  out. 

•  Special  Inputs:  No  comments. 

Tendency  to  overcontrol,  no  PIO. 

Had  to  work  at  it. 


•  PIO  Tendency 
AIRSPEEO  CONTROL: 


PERFORMANCE : 

t  Approach  Tasks: 

ILS:  Has  not  outstanding. 

Visual  (Sidestep):  Approaches  showed  a  slight  tendency  to  overcontrol. 

•  Landing  tasks:  Not  dor.e. 

•  Differences:  No  comments. 

WIND  AND  TURBULENCE:  Not  a  factor. 


SUH1ARY  COMMENTS:  Little  confused  about  this  one,  lot  of  distractions  dur 

ing  evaluations.  Not  a  solid  aircraft  but  not  bad,  un* 
certain  of  rating. 
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FEEL: 

•  Forces:  Okay. 

•  Displacement:  A  little  large. 

e  Sensitivity:  Okay,  detected  a  little  lag. 

PITCH  ATTITUDE  RESPONSE: 

e  Initial:  Initial  response  was  slow, 

e  Predictability:  Good. 

e  Special  Inputs:  No  special  pilot  techniques  (perhaps  a  little  lead.) 

e  PIO  Tendency:  No  PIO  tendency. 

AIRSPEEO  CONTROL:  Fairly  good. 

PERFORMANCE . 

e  Approach  Tasks: 

ILS:  Was  okay  (bad  pilot  performance.) 

Visual  (Sidestep):  Visual  approaches  fine, 
e  Landing  tasks:  Not  done, 

e  Differences:  Not  applicable. 

WIND  AND  TURBULENCE:  Smooth  air. 


SUMMARY  COMMENTS: 


Minor  problem  was  small  pitch  lag  and  a  little  more 
motion  than  desired.  Ceneraly  a  pretty  good  airplane. 


GAE/AA/80D-3 


FEEL: 

•  Forces;  Fairly  high  but  wouldn't  want  lighter. 

•  Displacement:  No  comments. 

•  Sensitivity;  No  comments 

PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Lagged. 

•  Predictability:  Wasn't  very  good  but  not  too  bad  either. 

•  Special  Inputs:  Necessary  to  try  to  anticipate,  reduce  aggressiveness. 

•  PIO  Tendency:  Tendency  to  overcontrol. 

AIRSPEED  CONTROL:  More  attention  required  because  of  lags  in  pitch  but 

never  bad. 

PERFORMANCE : 

•  Approach  Tasks: 

ILS:  Able  to  do  it. 

Visual  (Sidestep):  Approaches  went  fairly  well. 

•  Landing  tasks:  Not  performed. 

•  Differences:  Not  applicable. 

WIND  AND  TURBULENCE:  Not  a  factor. 

No  comments. 


SUWARY  COMMENTS: 


GAE/AA/80D-3 


1-3 

CONFIGURATION 

CJSP  ■  1.0  Jsp  *  0.74 

—  /0.2S  /  -- 

'  '  ^  1 

FLIGHT/PILOi 

1885/A 

PILOT  RATING: 

UINO/X-UIND: 

OVERALL  --  APPROACH  6  PIO  2 

0/04  TURBULENCE:  Light 

SP  6 

.  =  0.26  i 

6es 

FEEL: 

•  Forces:  High  initial  forces,  comfortable  in  steady  state. 


•  Displacement:  Normal. 


•  Sensitivity:  No  comments. 


PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Initial  response  sluggish. 


•  Predictability:  Final  response  not  predictable;  tendency  to  overshoot  a 

couple  of  tines. 

e  Special  Inputs:  Had  to  fly  with  snail  corrections  otherwise  easy  to 

overcontrol . 

e  PIG  Tendency:  No  comments. 


AIRSPEED  CONTROL:  Affected  by  control  technique  required  by  pitch. 


PERFORMANCE : 

e  Approach  Tasks: 

ILS:  Fairly  good,  had  to  pay  attention  to  pitch. 


Visual  (Sidestep):  Easy  to  ovcrcontrol  on  final  approach. 


•  Landing  tasks:  Nat  performed. 


•  Differences:  Not  applicable. 


WIND  AND  TURBULENCE:  Not  a  factor. 


SUMMARY  COMMENTS:  Tendency  to  overcontrol  in  pitch.  Difficult  to  make 

latge  corrections  fast. 
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FEEL: 

•  Forces;  Light  initially,  then  heavy;  selected  on  heavy  side. 

•  Displacement:  Felt  large,  spongy. 

•  Sensitivity:  Desirable  level. 

PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Very  delayed. 

•  Predictability:  Poor. 

e  Special  Inputs:  Keep  your  gain  down  and  fly  slowly  like  a  big  airplane. 

•  PIO  Tendency:  Definite  tendency,  tended  to  oscillate  2  or  3  cycles. 

AIRSPEED  CONTROL:  Bothersome,  poor,  objectionable. 

PERFORMANCE: 

e  Approach  Tasks: 

ILS:  Not  too  bad,  but  heavily  loaded. 

Visual  (Sidestep):  Felt  behind  the  airplane  on  sidesteps;  approaches 
better  than  on  IIS. 

•  Landing  tasks:  -  Not  done. 

•  Differences:  Not  applicable. 

U 1  NO  AND  TURBULENCE:  Moderate  corsswind  and  turbulence. 

SUMMARY  COMMENTS:  Hard  to  fly.  Very  delayed  initial  response  made  for  a 

high  workload,  particularly  on  the  ILS. 


GAE/AA/80D-3 


FEEL: 

•  Forces:  Forces  fairly  light,  greater  tendency  to  overcontrol  with 

heavier  forces. 

•  Displacement:  No  comments. 


•  Sensitivity:  No  comments. 


PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Initial  response:  some  lag  but  not  significant. 


•  Predictability:  Final  response  predictability  was  poor. 

•  Special  Inputs:  No  comments. 


•  PIO  Tendency:  No  comments. 

AIRSPEEO  CONTROL:  Reasonably  good. 


PERFORMANCE: 

•  Approach  Tasks: 

ILS:  Could  correct  large  errors,  some  tendency  to  PIO  on 

breakout. 

Visual  (Sidestep):  PIO  on  side  step  maneuver.  Tendency  to  PIO  on 
go-around. 

•  Landing  tasks:  Not  done. 

•  Differences:  Not  applicable,  but  worst  in  later  stage  of  approach. 


WIND  AND  TURBULENCE:  Crosswind  requires  extra  attention. 

SUMMARY  COMMENTS:  Difficult  to  acquire  an  attitude  quickly  and  predictably. 

Tendency  to  PIO. 
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FEEL: 

•  Forces:  Comfortable,  perhaps  a  bit  light. 

•  Displacement:  No  comments. 

•  Sensitivity:  No  comments. 

PITCH  ATTITUDE  RESPONSE: 

•  Initial:  Not  matched  with  response  after  input. 

•  Predictability:  Poor;  aircraft  bounces  in  turbulence  plus  a  tendency  to 

P10. 

•  Special  Inputs:  Had  to  spend  a  lot  of  time  on  attitude. 

•  PIO  Tendency:  Yes,  had  to  work  hard. 

AIRSPEED  CONTROL:  Not  bad  but  worked  hard. 

PERFORMANCE : 

•  Approach  Tasks: 

ILS:  Reasonable  but  hard  work. 

Visual  (Sidestep):  In  a  constant  oscillation  which  was  largely  pilot 
induced. 

•  Landing  tasks:  Not  done. 

•  Differences:  Definite  difference  between  instrument  .  isual 

approaches.  More  TIO  tendency  r  ...ui 

WIND  AND  TURBULENCE:  Aircraft  responded  which  compoun'.vd  difficulties. 


$ 


SUMMARY  COMMENTS: 


Tendency  to  bobble  and  to  PIO  were  major  problems. 
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Vita 

Randall  M.  Enright  was  born  on  18  May  1951  at  Scott  AFB,  Illinois. 

He  graduated  from  high  school  in  Del  Valle,  Texas  in  1969  and  attended 
the  United  States  Air  Force  Academy  from  which  he  received  the  degree  of 
Bachelor  of  Aeronautical  Engineering  in  June  1973.  He  completed  navigator 
training  and  received  his  wings  in  April  197^,  and  he  completed  electronic 
warfare  officer  training  in  October  197**.  He  served  as  a  B-52C-  Electronic 
Warfare  Officer  instructor /evaluator  with  the  51st  Bombardment  Squadron, 
Seymour  Johnson  AFB,  North  Carolina  until  May  1979-  During  that  time  he 
received  the  degree  of  Master  of  Business  Management  from  Central  Michigan 
University.  He  entered  the  School  of  Engineering,  Air  Force  Institute 
of  Technology,  in  June  1979-  He  is  presently  assigned  as  the  Lead  Stability 
and  Control  Engineer  for  the  Hypervelocity  Missile  Program,  Air  Force 
Armament  Laboratory,  Eglin  AFB,  Florida. 


Permanent  address:  1*721  Groveton  Way 

St.  Louis,  Mo.  63128 
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